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ABSTRACT 


The  use  of  finite  element  modeling  (FEM)  in  design  has  expanded  as 
computers  have  become  more  capable.  Despite  these  advancements,  the 
construction  of  physical  prototypes  remains  an  essential  aspect  of  design  and 
testing.  FEM  limitations  include  the  inability  to  accurately  account  for  joints, 
damping,  and  geometric  complexities.  Due  to  the  reality  gap  between  a  FEM  and 
the  prototype,  there  may  be  design  deficiencies  that  cannot  be  identified  until  the 
prototype  is  tested.  Using  eigenvalue  sensitivities,  enhanced  by  artificial 
boundary  conditions  (ABC),  the  gap  between  simulation  and  reality  can  be 
closed  via  FEM  updating.  With  an  updated  FEM,  the  same  eigenvalue 
sensitivities  can  be  utilized  to  detect  damage  in  structural  systems  in  use. 
Damage  that  produces  differences  in  natural  frequencies  between  the  structure 
and  its  FEM  can  be  related  to  the  loss  in  flexural  rigidity,  as  it  is  usually  assumed 
that  mass  modeling  is  correct.  This  indicator  allows  adjustment  of  a  FEM  to 
match  a  prototype  or  to  detect  damage  in  a  potentially  compromised  structure  via 
comparison  to  an  updated  FEM.  Based  on  simulation,  a  combination  of  multiple 
pin  and  spring  ABCs  is  optimal  for  producing  an  ideal  sensitivity  matrix,  and  thus, 
ideal  damage  detection  capability.  Flowever,  in  the  experimental  realm,  the 
synthesis  transformation  used  to  apply  ABCs  to  the  measured  frequency 
response  functions  can  distort  the  frequency  response  function  peaks,  leading  to 
error.  A  compromise  of  a  single  pin  ABC  permits  both  effective  model  updating 
and  damage  detection. 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  research  presented  here  centers  on  the  usage  of  the  finite  element 
method  (FEM).  Engineers  have  always  faced  problems  dependent  on  time  and 
space  throughout  the  design  process.  These  calculations  require  the  use  of 
partial  differential  equations  (PDEs),  which  may  be  unsolvable  via  analytic 
means.  Enter  the  FEM.  As  described  by  computer  modeling  company,  COMSOL, 
this  alternative  approach  involves  discretizing  the  PDE  to  create  an 
approximation  that  can  be  solved  using  numerical  methods  [1].  Approximations 
can  be  generated  through  the  linear  combination  of  base  functions  and  the 
application  of  boundary  conditions.  Figure  1  illustrates  this  process.  FEM  can  be 
applied  to  a  variety  of  scenarios,  such  as  convection-diffusion  problems,  Navier- 
Stokes  equations,  and  structural  dynamics  [1]. 


U 


i  i 


X 


Figure  1.  Discretized  Function  Approximation.  Source:  [1].  The  initial 
function  (blue)  is  seen  alongside  its  approximation  (red  dashes), 
which  is  a  linear  combination  of  basis  functions  (black). 

As  structural  dynamics  is  the  subject  of  concern  for  this  thesis,  it  is 
important  to  specify  that  the  linear  theory  of  elasticity  is  utilized.  This  permits  the 
linear  association  between  stress  and  strain.  It  is,  therefore,  essential  that  any 
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testing  does  not  result  in  exceeding  the  proportional  limit  of  the  test  sample’s 
stress-strain  curve.  Figure  2  illustrates  the  linear  stress-strain  relationship  below 
this  limit.  Obviously,  any  plastic  deformation  is  unacceptable.  With  these  criteria 
in  mind,  Euler-Bernoulli  beam  theory  may  be  properly  applied. 


"  (6) 


Figure  2.  Sample  Stress-Strain  Curve.  Source:  [2],  Note  the 
linear  stress-strain  relationship  preceding  the  proportional  limit. 


With  the  aid  of  Euler-Bernoulli  beam  theory,  the  beam  bending  problem  is 
simplified.  According  to  O.  A.  Bauchau  and  J.  I.  Craig,  this  theory  is  based  on 
three  fundamental  assumptions  [3].  First,  the  beam  cross  section  remains 
undeformed  in  its  own  plane.  Second,  the  cross  section  of  a  beam  remains  plane 
after  deformation.  Third,  the  cross  section  remains  normal  to  the  deformed 
neutral  axis  [3].  This  allows  for  relative  ease  in  calculating  the  response  of  a 
beam  when  subjected  to  various  loads. 

With  Euler-Bernoulli  beam  theory  at  one’s  disposal,  a  FEM  can  be 
generated  that  adheres  to  its  principles,  thus  permitting  the  simulation  of  load 
application  and  resultant  responses.  To  expedite  computation  time,  Euler- 
Bernoulli  beam  theory  is  applied  using  simple  one-dimensional  line  model  to 
represent  the  beam  in  question.  Due  to  the  beam’s  uncomplicated  geometry,  the 
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use  of  such  a  basic  FEM  is  acceptable.  An  additional  consideration  is  the  level  of 
discretization  in  the  simulated  beam.  As  seen  in  Figure  1,  the  function  is  divided 
into  discreet  components,  each  of  which  is  approximated  using  a  linear 
combination  of  basis  functions.  Likewise,  the  beam  FEM  can  be  discretized  into  a 
certain  number  of  elements.  A  high  element  count  can  provide  approximations 
that  are  more  accurate.  However,  as  the  number  of  elements  increases,  so  does 
the  computation  time.  This  processing  time  is  a  major  limiting  factor  of  FEM 
complexity.  As  with  function  u  in  Figure  1 ,  the  size  of  the  elements  can  be  varied, 
so  relatively  uniform  geometry  is  covered  by  larger  elements,  while  sharper 
changes  in  geometry  are  accounted  for  with  smaller  elements. 

While  increasing  the  element  count  can  improve  the  accuracy  of  the  FEM, 
it  does  not  entirely  eliminate  inconsistencies  between  a  simulated  item  and  its 
physical  counterpart.  There  remain  unique  complexities  in  the  real  world  that  are 
not  initially  accounted  for  by  FEM.  For  example,  structural  joints  and  damping  are 
phenomena  that  cannot  be  accurately  calculated  without  experimental  input  and 
will,  therefore,  be  lacking  in  a  first  iteration  FEM.  Additionally,  physical 
irregularities  such  as  deviations  from  material  homogeny,  variations  in 
production,  etc.,  may  exist.  These  differences  between  simulation  and  reality  can 
be  revealed  by  generating  and  comparing  frequency  response  functions  (FRF) 
for  the  FEM  and  its  physical  associate. 

Of  particular  interest  is  the  frequency  at  which  function  peaks  exist.  The 
peaks  of  an  FRF  generated  from  simulated  data  correspond  to  natural 
frequencies.  FRFs  obtained  via  experimental  testing  can  also  identify  natural 
frequencies,  once  damping  is  accounted  for.  By  comparing  these  two  functions,  a 
matrix  of  frequency  differences  can  be  assembled.  This,  in  turn,  can  be  applied 
to  eigenvalue  sensitivities  to  determine  the  structural  discrepancy.  Figure  3 
exhibits  a  typical  FRF. 
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Figure  3.  Example  FRF  Generated  from  Experimental  Data 

B.  LITERATURE  REVIEW 

The  concept  of  using  artificial  boundary  conditions  for  model  updating  and 
damage  detection  has  been  explored  by  J.  H.  Gordis  since  the  1990s.  His  article 
“Artificial  Boundary  Conditions  for  Model  Updating  and  Damage  Detection”  in 
Mechanical  Systems  and  Signal  Processing  describes  the  approach  [4],  This 
published  work  states  that  FEM  improvement  in  required  to  establish  an 
acceptable  level  of  confidence  in  the  FEM’s  ability  to  predict  structure  response. 
According  to  Gordis,  “This  inaccuracy  of  the  FE  model  is  reduced  (the  model  is 
‘improved’)  by  the  adjustment  of  selected  physical  and  material  parameters 
which  define  the  model.  These  parameters  can  include  dimensional  properties  of 
structural  elements,  moduli  of  elasticity,  and  densities,  for  example”  [4],  An 
inherent  complication  with  this  approach  is  the  fact  that  modal  tests  yield  a 
number  of  modal  parameters  that  is  usually  less  than  the  number  of  parameters 
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that  define  the  FEM.  This  produces  an  underdetermined  problem,  as  there  is  not 
enough  usable  data  to  reliable  adjust  (update)  the  FEM  [4], 

In  Reference  [4],  Gordis  notes  the  work  of  S.  Li,  S.  Shelley,  and  D.  Brown, 
who  introduced  a  method  of  expanding  the  measured  modal  parameter  database 
through  the  use  of  perturbed  boundary  conditions  (PBCs).  This  approach 
requires  physical  modifications  to  the  test  structure.  This  can  include  the  addition 
of  mass  at  selected  points.  With  every  different  modification,  an  additional  modal 
test  must  be  conducted  [4],  While  this  method  does  provide  the  requisite  data,  it 
is  relatively  time  consuming  and  costly  due  to  the  implementation  of  multiple 
structure  modifications  and  the  performance  of  numerous  modal  tests. 

Gordis  presents  an  alternative  solution.  Through  the  use  of  artificial 
boundary  conditions  (ABCs)  vice  PBCs,  physical  modifications  can  be  avoided. 
Imposing  boundary  conditions,  such  as  the  mass  additions,  can  be  easily  done  in 
the  simulated  realm.  Though  these  ABCs  produce  multiple  FEMs,  corresponding 
measured  mode  frequencies  can  be  obtained  from  a  single  modal  test. 
Additionally,  more  complex  modifications,  such  as  the  application  of  pins,  can  be 
applied,  thus  creating  an  even  larger  database  of  modal  parameters  [4],  As 
stated  in  a  subsequent  article  by  J.  H.  Gordis  and  K.  Papagiannakis,  for  each 
pin-based  ABC,  an  omitted  coordinate  system  (OCS)  is  defined.  This  system 
corresponds  to  a  structure  with  additional  restraints  placed  at  some  combination 
of  locations  [5]. 

Gordis  reveals  that  without  the  application  of  boundary  conditions,  the 
parameter  database  remains  limited,  and  the  eigenvalue  sensitivity  matrix 
(assembled  from  the  parameter  database)  risks  containing  data  sets  that  are 
linearly  dependent  (or  near-linearly  dependent)  to  one  another  [4],  This  linear 
dependence  results  in  a  poorly  conditioned  matrix.  In  other  words,  it  is  a  matrix 
that  is  rank  deficient.  The  unfortunate  result  is  a  rank  deficient  matrix  that  can  fail 
to  produce  the  required  adjustments  for  model  updating  or  detect  damage  in 
every  element. 
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Even  with  a  large  database  of  modal  parameters  generated  from  both 
baseline  FEMs  and  those  with  ABCs  applied,  a  well-conditioned  sensitivity  matrix 
is  not  guaranteed.  An  ideal  combination  of  parameters  must  be  selected  for 
sensitivity  matrix  assembly.  Gordis  and  Papagiannakis  offer  an  approach  to 
sensitivity  matrix  construction  using  QR  decomposition  with  column  pivoting  [5]. 
In  their  approach,  all  pin-based  ABC  location  combinations  are  applied  to  a 
simple  beam  FEM.  Errors  were  intentionally  introduced  to  the  FEM,  and  each 
parameter  set  was  utilized  in  an  attempt  to  localize  these  errors.  In  each 
scenario,  relative  magnitude  error  (RME)  and  maximum  false  alarm  (MFA)  were 
calculated.  When  these  criteria  reached  near-zero  values,  the  sensitivity  row  set 
was  saved.  The  specific  columns  selected  for  this  ideal  set  are  identified  by 
utilizing  the  permutation  matrix  of  the  QR  decomposition  [5]. 

J.  Shin  and  J.  FI.  Gordis  offer  a  simplified,  alternative  approach  to 
parameter  selection  called  orthogonal  projection.  This  method  is  focused  on 
addressing  the  condition  number  of  the  sensitivity  matrix.  The  condition  number 
as  a  direct  indication  of  the  level  of  linear  dependence  between  data  sets  within 
the  matrix.  The  lower  the  condition  number,  the  less  linear  dependence  exists. 
With  linear  dependence  minimized,  model  update  and  damage  detection 
accuracy  is  maximized  [6].  This  method  can  be  universally  applied  to  any  ABC 
type,  location,  or  quantity  without  modification.  Due  to  its  simplicity  and  ease  of 
implementation,  this  approach  is  preferred  for  use  in  this  thesis. 

C.  SCOPE  OF  THESIS 

This  thesis  delves  into  the  variety  of  ABCs  that  can  be  used  to  populate  a 
candidate  matrix  for  data  set  selection  and  eigenvalue  sensitivity  optimization. 
Different  types,  numbers,  and  locations  of  ABCs  are  explored  in  an  effort  to 
provide  a  substantial  library  of  parameter  and  modal  options. 

Additionally,  different  styles  of  FEMs  are  compared  with  regard  to 
accuracy  and  computation  time.  The  structure  that  is  modeled  is  a  rectangular 
cross-section  free-free  beam.  Since  force  application  and  resultant  translation  is 
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uniaxial,  a  simple  two-DOF  per  node  FEM  is  utilized.  Flowever,  different  element 
sizes  and  numbers  can  be  tried.  These  various  models  are  used  to  identify  a 
candidate  “damaged  condition”  for  the  test  beam. 

Finally,  experimental  validation  of  the  eigenvalue  sensitivity  concept  is 
investigated.  Associated  complications  with  working  in  physical  realm  will  be 
discussed  along  with  devised  solutions.  Recommendations  for  achieving  the  best 
results  are  presented. 

D.  SOFTWARE  NOTES 

FEM  construction,  model  updating,  and  simulated  damage  detection  is 
performed  using  MATLAB  2017a  software  [7],  Beam  impact  testing  and  data 
collection  is  conducted  with  Pulse  Reflex  software  by  Bruel  &  Kjaer  [8]. 
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II.  THEORY 


A.  EULER-BERNOULLI  BEAM  THEORY 

Since  all  FEMs  constructed  for  this  thesis  follow  Euler-Bernoulli  beam 
theory,  it  is  important  that  the  associated  equations  are  discussed.  This  begins 
with  considering  the  forces  and  moments  acting  on  a  beam.  W.  T.  Thomson  and 
M.  D.  Dahleh  state  that  the  sum  of  forces  is  the  /-direction  is  [9]: 

dV-p(x)dx  =  0  (1) 

where  V  is  the  shear  force  and  p(x)  is  the  loading  per  unit  length  of  the  beam. 
The  sum  of  moments  about  any  point  is  [9]: 

dM  -  Vdx  -  j  p(x)(dx)2  =  0  (2) 

where  M  is  the  bending  moment  [9].  Figure  4  illustrates  these  variables. 


These  equations  provide  the  following  relationships  [9]: 
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dV 

dx 


P(x) 


dM 

dx 


=  V 


Based  on  these  relationships,  the  following  equation  is  obtained  [9]: 


drM 

dx2 


dV_ 

dx 


p(x) 


The  bending  moment,  M,  is  related  to  the  flexure  equation  [9]: 


(3) 

(4) 

(5) 


M  =  El 


Substituting  /Win  Equation  (5)  produces  [9]: 


(6) 


d2 

dx2 


f 

El 

v 


1 2  'N 

a  y 
dx2  j 


P(x) 


(7) 


According  to  William  Thomson  and  Marie  Dahleh’s  Theory  of  Vibration 
with  Applications,  “beam  vibrating  about  its  static  equilibrium  position  under  its 
own  weight,  the  load  per  unit  length  is  equal  to  the  inertia  load  due  to  its  mass 
and  acceleration.  Because  the  inertia  force  is  in  the  same  direction  as  p(x)... we 
have,  by  assuming  harmonic  motion”  [9]: 


d2 

dx2 


f 

El 

V 


1 2  \ 

d  y 

dx2  J 


-  pro2  y 


0 


(8) 


where  E  is  the  modulus  of  elasticity,  /  is  the  second  moment  of  inertia,  p  is  mass 
density,  and  to  is  the  natural  frequency  [9]. 


B.  FINITE  ELEMENT  MODELING 

With  the  principles  of  Euler-Bernoulli  beam  theory  specified,  the  FEM  can 
be  constructed  accordingly.  For  a  simple,  one-dimensional  model,  each  node 
features  two  degrees  of  freedom  (DOF),  one  translational,  the  other  rotational. 
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According  to  R.  R.  Craig,  Jr.,  for  a  uniform  beam  element,  transverse 
displacement  is  defined  as  [10]: 


V(x,t)  =  YJV/i(t)vi(t)  (9) 

*=  1 

where  V  is  the  transverse  displacement  at  a  particular  point  along  the  length  of 
the  beam  element(x)  at  a  particular  moment  in  time  ( t ).  v  is  the  transverse  and 
angular  displacements  at  each  end  of  the  beam  element,  as  seen  in  Figure  5 
[10]- 


Figure  5.  Transverse  Displacement  of  a  Uniform  Beam  Element. 

Source:  [10]. 


V  is  the  shape  functions  (i.e.,  basis  functions)  that  are  linearly  combined  to 
replicate  the  transverse  displacement  function.  These  shape  functions  must 
adhere  to  the  following  boundary  conditions  [10]: 

Wx  (0)  =  1.  K(0)  =  Vi  (L)  =  vW  =  0 

^2(0) =1.  ^2(0)  =  ^2(L)  =  ^;(L)=o 

y/3(L)  =  \,  y/3(0)  =  ,/3(0)  =  ,//'?>(L)  =  0 
VW  =  1>  ^4(0)  =  ^4(0)  =  ^4(L)=0 
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Figure  6  illustrates  these  shape  functions  along  an  element  with  length  L 
and  their  respective  boundary  conditions. 


Figure  6.  Shape  Functions  Showing  Boundary  Conditions. 

Source:  [10]. 


The  general  solution  to  Equation  (9)  is  a  cubic  polynomial  [10]: 

V  (x)  =  q  +  c2 1 


M 

M 

2 

|  +  c,  | 

U  y 

+  C4 

(11) 


Substituting  the  boundary  conditions  into  Equation  (11)  produces  the  four  shape 
functions  [10]: 
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V\  =  1  —  3 1  y  +2  [j 
L  J  V  L 


¥l-x  2L\  J  +  *<(_) 


(12) 


f  x] 

2 

I 

-2 

I 

f  xY 

-  +L\ 

- 

UJ 

L 

Using  the  principles  of  Bernoulli-Euler  beams  along  with  these  shape 
functions,  stiffness  and  mass  coefficients  can  be  calculated  [10]: 


L 

kij  =  J  ElyA/Yclx 
0 

L 

m(j  =  |  pAy/xj/ylx 
0 


(13) 

(14) 


By  substituting  the  shape  functions  from  Equation  (12)  into  Equations  (13) 
and  (14),  the  element  stiffness  and  mass  matrices  can  be  obtained  [10]: 


K 


elem 


M 


elem 


(  El\ 


'  pAL 
v^20" 


12 

6  L 

-12 

6  L 

4  L2 

-6  L 

21} 

12 

-6  L 

symm. 

41} 

"  156 

22  L 

54 

-13  L 

4  L2 

13  L 

-31} 

156 

-22  L 

symm. 

41} 

(15) 


(16) 


It  is  important  to  note  that  these  matrices  only  account  for  a  single  element 
in  a  model  composing  of  N  elements,  N+1  nodes,  and  2*(N+1)  total  DOF.  All 
elements  must  be  adjoined  to  assemble  the  model  in  its  entirety  and  to  account 
for  the  fact  that  interior  elements  have  common  nodes  (and  their  respective  DOF) 
with  adjacent  elements.  This  is  done  by  overlapping  element  stiffness  and  mass 
matrices  at  DOFs  in  common  as  shown  on  the  following  page: 
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With  stiffness  and  mass  matrices  accounting  for  the  entire  FEM  (i.e., 
global  matrices),  the  equation  of  motion  can  be  applied  as  presented  by  W.  T. 
Thomson  and  M.D.  Dahleh  [9]: 

[M\{x}+[C\{x}+[K]{x}  =  {f(t)}  (19) 

where  [M\  is  the  global  (system)  mass  matrix,  [C]  is  the  global  damping  matrix, 
and  [K]  is  the  global  stiffness  matrix.  The  u  vectors,  from  left  to  right,  are 
transverse  acceleration,  transverse  velocity,  and  transverse  position.  Vector  f 
contains  the  external  forces  and  moments.  This  equation  of  motion  is  generally  a 
set  of  M  coupled  equations.  Fortunately,  through  the  use  of  modal 
transformation,  these  equations  can  be  decoupled  by  premultiplying  Equation 
(19)  by  the  modal  matrix  P  (the  modal  matrix  is  determined  by  the  solution  of  the 
homogeneous  undamped  equation).  Letting  {x}=[P\{y},  the  following  equation  is 
obtained  [9]: 

[pf  [M]lP]{y}  +  lPf  [C][P]W  +  [pf  [K][P]{y}  =  [Pf  {fit)}  (20) 

The  mass  and  stiffness  terms  of  Equation  (20)  produce  diagonal  matrices. 
However,  the  damping  term  is  generally  not  diagonal.  In  order  to  make  this  term 
diagonal,  proportional  damping  must  be  assumed.  These  actions  permit  the 
uncoupling  of  Equation  (19),  with  its  /th  equation  taking  the  following  form  [9]: 

y,  +  2 4m.  y, + y,  =  r,f,  (0  (2 1 ) 

where  £  is  the  damping  ratio,  defined  by  R.  R.  Craig,  Jr.,  is  [1 0]: 

€  =  —  (22) 

Ccr 

where  c  is  the  damping  coefficient  (obtained  experimentally)  and  ccr  is  the  critical 
damping  coefficient  [10]: 


ccr  -  2  mcoi  =  —  =  2 \fian 


co . 
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(23) 


6U/  is  natural  frequency,  which  is  defined  as  [10]: 


CO,  = 


(24) 


and  r,  is  the  mode  participation  factor,  which  is  defined  by  Thomson  and  Dahleh 
as  [9]: 


I ]Mxj)[p](xj) 

r,-  =  J -  (25) 

XmM2(-r) 

j 

Note  the  lack  of  an  element  and  global  damping  matrix  in  the  previous 
calculations.  The  experimental  rig  supporting  the  free-free  beam  in  question 
results  in  insignificant  damping.  Therefore,  [C]  can  be  omitted  from  Equation  (19) 
and  will  not  be  considered  when  simulating  dynamic  response.  Additionally,  the 
beam  will  not  be  subjected  to  external  loading.  Therefore,  Equation  (19)  can  be 
reduced  to: 


[M]{x}+[K]{x}  =  0  (26) 

However,  it  is  important  to  note  that  experimental  data  analysis  involves 
using  an  algorithm  that  calculates  damping  ratio  based  on  FRF  shape.  The 
relationship  between  damped  and  undamped  natural  frequency  is: 

o>d=ojJ\-f-  (27) 

where  cod  is  the  damped  natural  frequency,  ujn  is  the  undamped  natural 
frequency,  and  £is  the  damping  ratio. 

C.  EIGENVALUE  SENSITIVITY 

In  order  to  allow  for  model  updating  and  damage  detection,  a  relationship 
between  differences  in  natural  frequencies  and  differences  in  physical 
parameters  must  be  defined.  According  to  J.  H.  Gordis,  the  sensitivity  is  a  rate  of 
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change  of  some  eigenvalue,  A,  of  mode  /,  with  respect  to  the  change  of  some 
physical  parameter,  p  as  shown  here  [11]: 


(28) 


Therefore,  changes  in  an  eigenvalue  can  be  predicted  given  a  change  in  a 
physical  parameter  [11]: 


%  =  * ,+tVAPj)  (29) 

where  A) is  the  altered  eigenvalue  that  is  the  result  of  a  parameter  modification, 
Apj.  By  constructing  a  sensitivity  matrix  [S\  of  S,y  elements,  each  row  representing 
a  mode  and  each  column  representing  a  parameter,  the  following  equation  can 
be  formed  [11]: 


{A2}  =  [5]{Ap}  (30) 

Equation  (30)  is  what  allows  for  model  updating  by  comparing  simulated 
and  experimental  obtained  natural  frequencies.  It  is  also  what  provides  indication 
of  damage  in  the  form  of  Ap,  which  is  generated  from  a  variance  in  natural 
frequencies  from  an  updated  model  and  a  damaged  physical  beam  [11], 

To  generate  the  sensitivity  matrix,  one  must  calculate  each  element,  S,y.  To 
do  so,  the  eigensystem  must  first  be  defined  [11]: 

[K-XiM]{^}  =  {  0}  (31) 

where  0®  is  the  eigenvector  for  mode  /.  Both  [K]  and  [M\  are  functions  of  certain 
physical  parameters.  In  this  work,  the  perturbation  {Ap)  will  be  defined  only  as 
changes  in  flexural  rigidity  {E*l).  Therefore,  only  [K]  will  be  affected  by  these 
perturbations.  Taking  the  derivative  of  Equation  (31)  yields  [11]: 

a K  .  dM  dX, 

dPj  dp,  dP  j 


fS  o 


(32) 


17 


Pre-multiplying  by  f0^rand  transposing  Equation  (32)  produces  [11]: 


dK 

dM 

t: 

8A  " 

- -M 

+  < 

dPj 

SPj 

dPj 

Sp 

K  J  J 

0 


(33) 


Invoking  Equation  (29)  and  assuming  mass-normalized  shapes,  {(p(i>}T[M]  {(p(i>}=1, 
the  equation  can  be  manipulated  to  yield  [11]: 


dA,. 

Spj 


--A  — 
dPj  '  dPj 


0 


(34) 


It  is  essential  to  note  that  these  derivatives  are  calculated  by  using  a 
forward  difference  approximation  based  on  the  linear  terms  of  the  Taylor  series 
[11].  Therefore,  if  8A,  from  Equation  (28)  exceeds  the  linearity  of  the  eigenvalue 
sensitivity,  the  equation  will  not  be  usable. 


D.  APPLICATION  OF  ARTIFICIAL  BOUNDARY  CONDITIONS  IN  THE 

FINITE  ELEMENT  MODEL 

Due  to  linear  dependence  between  the  rows  of  the  eigenvalue  sensitivity 
matrix,  the  rank  of  the  matrix  is  less  than  the  number  of  parameters.  This  can 
lead  to  a  failure  to  detect  damage  in  certain  elements.  In  order  to  bolster  the 
capability  of  the  sensitivity  matrix,  ABCs  are  applied  to  allow  the  construction  of  a 
full-rank  sensitivity  matrix.  However,  a  full  rank  sensitivity  matrix  is  only  a 
minimum  standard  of  acceptability.  J.  Shin  and  J.  H.  Gordis  state  that  effective 
model  updating  and  accurate  damage  detection  is  also  dependent  on  the 
matrix’s  condition  number  [6].  This  number  can  be  driven  toward  unity  (perfect 
linear  independence)  with  the  judicious  use  of  ABCs.  The  types  of  ABCs 
examined  here  include  artificial  pins,  point  masses,  and  springs. 

Pins  eliminate  translation  while  preserving  rotation  at  their  location. 
Therefore,  an  artificial  pin  will  eliminate  the  translational  DOF  at  the  FEM  element 
node  at  which  it  is  applied.  This  effect  is  represented  in  the  FEM  by  eliminating 
the  row  and  column  associated  with  the  pinned  DOF  in  the  global  mass  and 
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stiffness  matrices.  For  springs  to  ground,  the  diagonal  stiffness  matrix  element 
associated  with  the  modified  DOF  is  increased  by  the  spring  constant.  Point 
masses  are  applied  in  a  similar  fashion,  except  it  is  the  diagonal  mass  matrix 
element  that  is  adjusted.  Like  the  pins,  springs  and  point  masses  only  affect  the 
translational  DOFs;  rotational  DOFs  remain  unaltered.  Multiple  ABCs  can  be 
applied  at  various  nodes.  A  key  consideration  is  that  the  ABC  locations  on  the 
FEM  must  correspond  to  the  accelerometers  placed  on  the  physical  structure. 
This  way,  the  same  ABC  configuration  can  be  applied  to  the  measured  data. 

Each  ABC  configuration  provides  its  own  sensitivity  matrix.  These 
matrices  can  be  stacked  into  one  all-encompassing  sensitivity  library  matrix  as 
shown  here: 


[s»>] 

[s'11] 


library 


E.  SELECTION  OF  ARTIFICIAL  BOUNDARY  CONDITIONS 


(35) 


J.  Shin  and  J.  H.  Gordis  state  that  with  an  eigenvalue  sensitivity  matrix 
library  established,  individual  rows  from  this  matrix  must  be  selected  to  form  a 
square  matrix  that  can  be  applied  in  Equation  (30)  [6]: 


K1} 

1 

1  i 

§ 

Co 

1 _ 1 

K1} 

CN 

V  •• 

>  = 

K] 

K'>} 

KT 

(M 


(36) 


19 


These  matrices  are  referred  to  as  “composites”  as  results  in  Equation  (36)  being 
rewritten  as  [6]: 


{AA}=[S]{A,,}  (37) 

where  AA  is  the  composite  vector  of  eigenvalue  differences,  and  S  is  the 

composite  sensitivity  matrix  [6].  Each  row  from  the  eigenvalue  sensitivity  matrix 
library  is  a  candidate  for  populating  the  composite  sensitivity  matrix.  These 
individual  rows  are  identified  as  sS,y  where  L  is  the  total  number  of  rows  in  the 
library,  /  is  the  row  index  and  j  is  the  ABC  configuration.  Normalized  rows  are 
defined  as  [6]: 


ft/ 


(38) 


The  main  objective  in  populating  the  composite  sensitivity  matrix  is  to 
append  rows  from  the  candidate  library  that  maximize  the  angle  between  these 

candidates  and  the  subspace  spanned  by  previously  selected  rows,  R(S  ). 

During  the  composite  matrix  assembly,  the  set  of  rows  selected  thus  far  is  S 

and  the  candidate  row  is  ft; .  The  angle  between  and  S  is  produced  by  [6]: 


{%} 

ft 

• 

where  [P\  is  the  projector  defined  as  [6]: 


(39) 


P  =  s(sTs)lST  (40) 

Figure  7  illustrates  the  vector  of  a  candidate  row  projecting  onto  the 
subspace  spanned  by  rows  already  selected. 
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Figure  7.  Candidate  Row  Vector  and  Selected  Row  Subspace. 

Source  [6]. 

In  addition  to  the  condition  number,  the  mode  dependency  index  (MDI) 
and  the  parameter  dependency  index  (PDI)  also  serve  as  indicators  of  linear 
dependency.  Specifically,  MDI  and  PDI  reflect  the  linear  independency  between 
any  two  rows  or  two  columns  in  the  sensitivity  matrix,  respectively.  MDI  and  PDI 
are  defined  as  [6]: 


MDI,.  = 


PD‘uj  = 


Mfr J 
MfcJ 


(41) 


(42) 


where  R  is  one  of  the  two  rows  of  the  composite  sensitivity  matrix  being 
examined  for  linear  dependence,  and  C  is  one  of  the  two  columns  [6]. 

F.  FREQUENCY  RESPONSE  FUNCTION  ANALYSIS 

With  the  composite  eigenvalue  sensitivity  matrix  and  eigenvalue  difference 
vector  populated,  the  associated  ABC  configurations  must  now  be  applied  to 
experimental  data.  Modal  testing  of  a  physical  beam  does  not  yield  global  mass 
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and  stiffness  matrices,  such  as  the  ones  generated  by  the  FEM.  Therefore,  a 
different  approach  is  required  to  apply  the  selected  ABC  configurations  to  the  test 
data. 

According  to  J.  H.  Gordis,  the  FRF  matrix,  [H\,  is  generated  by  the  modal 
test.  Ignoring  damping,  each  FRF  matrix  is  defined  as  [4]: 


m 

co;-Qr 


(43) 


where  is  0r,a  mass  normalized  mode  shape,  wris  the  rth  natural  frequency,  and 
Q  is  the  forcing  frequency.  These  FRFs  are  then  assigned  to  different  coordinate 
sets.  The  first  set,  the  connection  set  (CSET),  contains  FRFs  that  belong  to  a 
modified  DOF.  For  example,  the  placement  of  an  artificial  pin  will  result  in  the 
translational  DOF  of  the  pinned  node  being  assigned  to  the  CSET.  The 
remaining  DOFs  (not  altered  by  ABCs)  are  assigned  to  the  internal  set  (ISET). 
The  partitioned  FRF  matrices  take  the  following  form  [4]: 


H,, 

Hic 

H. 

H . 

(44) 


The  new  FRF  can  be  generated  by  applying  the  partitions  to  the  following 
formula  [4]: 


H*=Hii-HicH2Hci  (45) 

This  modified  FRF  matrix,  H*,  is  recalculated  for  every  frequency  in  the 
desired  bandwidth.  Taking  the  average  of  its  diagonals  will  yield  an  FRF  that  all 
of  the  modes  within  the  test  bandwidth.  Plotting  this  function  will  reveal  amplitude 
peaks  that  correspond  to  the  test  article’s  natural  frequencies.  These  frequencies 
can  then  be  compared  to  the  FEM  natural  frequencies,  thus  producing  the 
composite  vector  of  eigenvalue  differences.  This  permits  the  use  of  Equation  (37) 
to  solve  for  any  parameter  change  required  for  model  updating  or  for  detecting 
damage  [4], 
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III.  EXPERIMENTAL  APPROACH 


A.  TEST  ARTICLES 

The  test  articles  utilized  for  modal  testing  are  extruded  aluminum 
rectangular  beams.  Upon  receipt,  these  beams  were  drilled  and  tapped  to 
accommodate  accelerometer  mounting  screws.  The  first  beam  was  not  altered 
any  further  to  be  used  as  an  “intact”  structure.  The  second  was  reserved  for 
further  alterations  to  act  as  damage.  Table  1  lists  the  beam  properties  and 
Appendix  A  features  the  inspection  report. 


Table  1 .  Test  Beam  Properties 


Material 

Al  6061 -T651 1 

Mass,  kg  (Ibm) 

2.449  (5.4) 

Length,  cm  (in) 

182.88  (72) 

Width,  cm  (in) 

5.09  (2.004) 

Thickness,  cm  (in) 

0.9627  (0.379) 

Tap  Diameter,  cm  (in) 

0.4064  (0.16) 

Distance  from  Beam  Edge  to  Tap,  cm  (in) 

1.27  (0.5) 

Distance  Between  Taps,  cm  (in) 

12.8778  (5.07) 

Due  to  the  focus  on  ABCs,  the  beam  was  suspended  in  a  free-free 
condition.  This  was  accomplished  by  utilizing  two  elastic  cords  attached  to  an 
overhead,  each  about  270  cm  (106.3  in)  in  length.  The  resultant  pendulum 
frequency  is  0.54  Hz,  far  less  than  any  natural  frequency  and,  therefore,  does  not 
interfere  with  the  modal  testing.  Figure  8  shows  one  of  the  rectangular  beams 
and  Figure  9  shows  the  test  rig  with  the  suspended,  free-free  test  article. 


Figure  8.  Rectangular  Beam  With  Fifteen  Taps  for  Accelerometers 
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Figure  9.  Suspended  Beam 

B.  INSTRUMENTATION 

Each  of  the  fifteen  taps  allows  for  the  attachment  of  two  accelerometer 
10-32  mounting  studs  (shown  in  Figure  10).  One  is  used  for  attaching  the 
accelerometer,  while  the  other  is  used  as  an  impact  point  for  the  force  hammer. 


Figure  10.  Mounting  Studs  with  2.54  cm  (1  in)  Scale 
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The  accelerometers  are  PCB  Piezotronics  ICP  Model  336C04  (shown  in  Figure 
11).  Table  2  contains  the  dimensions  of  the  accelerometer  and  its  mount. 


Figure  1 1 .  Accelerometers  with  5.08  cm  (2  in)  Scale 


Table  2.  Accelerometer  and  Mount  Properties 


Accelerometer  Mass,  g  (oz) 

5.07  (0.179) 

Accelerometer  Diameter,  cm  (in) 

1.575  (0.62) 

Accelerometer  Height,  cm  (in) 

1.397  (0.55) 

Mounting  Screw  Mass,  g  (oz) 

0.72  (0.0254) 

Mounting  Screw  Diameter,  cm  (in) 

0.4191  (0.165) 

Mounting  Screw  Height,  cm  (in) 

0.7061  (0.278) 

The  impact  hammer  is  a  PCB  Piezotronics  Model  086B03  (shown  in  Figure  12). 
Before  conducting  the  modal  test,  each  accelerometer’s  calibration  was  verified 
using  a  PCB  Piezotronics  Model  194C06  handheld  shaker  (shown  in  Figure  13) 
in  conjunction  with  the  Pulse  Reflex  Software.  Appendix  B  features  the 
calibration  and  specification  sheets  for  each  of  these  instruments. 
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Figure  12.  Impact  Hammer 


Figure  1 3.  Handheld  Shaker  with  Accelerometer  Attached  at  Top 
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C.  MODAL  TESTING 

Modal  testing  can  be  conducted  in  multiple  fashions,  each  with  their  own 
unique  method  of  providing  beam  excitation:  step  sine,  random  burst,  single 
impact,  and  random  impact. 

The  step  sine  approach  utilizes  a  shaker  and  force  transducer  (load  cell) 
attached  to  a  mounting  screw.  The  shaker  must  be  suspended  from  a  stand  to 
maintain  the  rig’s  free-free  condition.  An  additional  consideration  is  the  length 
and  stiffness  of  the  attachment  rod  (“stinger”).  A  stinger  that  is  too  long  or  lacks 
the  adequate  stiffness  will  result  in  energy  dissipation  and  skew  results.  Beam 
excitation  is  conducted  at  a  specified  frequency  interval  across  the  entirety  of  the 
desired  bandwidth.  The  Pulse  Reflex  software  pauses  at  each  frequency  to 
ensure  that  all  transients  have  dissipated  and  that  beam  response  is  caused 
solely  by  excitation  at  the  dwell  frequency.  The  higher  the  frequency  resolution, 
the  longer  the  process  takes  to  cover  a  specified  bandwidth.  Due  to  the  minimum 
desired  frequency  resolution  of  0.3125  Hz  and  the  1000  Hz  bandwidth,  the  step 
sine  approach  is  extremely  time  consuming. 

The  burst  random  technique  also  employs  a  shaker  and  accompanying 
force  transducer.  As  the  name  suggests,  the  shaker  generates  a  random  time 
signal  that  contains  all  frequencies  within  the  bandwidth  of  interest  to  produce  the 
beam’s  dynamic  response.  Compared  to  step  sine,  this  method  is  significantly 
quicker  to  complete. 

The  single  impact  method  uses  a  force  hammer  to  impact  each  node 
where  an  accelerometer  is  mounted.  The  node  is  struck  once  with  the  hammer 
and  the  response  is  recorded  for  a  period  of  time.  The  higher  the  desired 
frequency  resolution,  the  longer  the  recording  time  required.  To  ensure 
consistent  results,  these  single  impacts  are  repeated  several  times  and  the 
resulting  frequency  response  functions  are  averaged  at  each  node. 

Random  impact  also  utilizes  the  impact  hammer.  Instead  of  using  a  single 
hammer  strike,  numerous  hammer  impacts  are  applied  over  a  period  of  time.  The 
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times  at  which  data  is  collected  by  the  Pulse  Reflex  software  are  randomized  in 
order  to  ensure  there  are  no  artificialities  produced  from  inadvertent  periodicity  in 
hammer  impacts. 

Out  of  the  four  methods,  the  single  and  random  impact  modal  tests  were 
deemed  most  useful  for  the  purposes  of  this  research.  The  shaker-based 
techniques  required  additional  setup  due  to  the  more  complex  rig.  Additionally, 
due  to  the  beam’s  low  mass,  the  attached  shaker  added  an  unacceptably  high 
mass  percentage  to  the  test  article,  which  could  potentially  corrupt  results.  As 
can  be  seen  in  Figure  14,  the  shaker-based  approaches  also  generated 
additional  noise  at  lower  frequencies.  The  hammer  methods,  on  the  other  hand, 
require  less  setup,  and  produce  a  cleaner  FRF  at  lower  frequencies.  In  addition, 
there  is  no  risk  of  mass  loading  of  the  test  article  by  the  shaker.  It  is  important  to 
note  that  the  single  impact  method,  while  effective  at  the  initially  utilized  0.3125 
Hz  frequency  resolution,  required  a  twelve  second  recording  time  for  a  resolution 
of  0.08  Hz.  Due  to  this  long  recording  time,  random  impact  was  utilized  for 
resolutions  of  0.08  Hz  and  higher. 
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Random  Impact 
Random  Burst 
Single  Impact 
Step  Sine 


Figure  14.  Experimental  FRF  Comparison.  Source:  [12]. 
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IV.  FEM  MESHING  AND  DAMAGE  SCENARIO  CANDIDACY 


The  first  issue  that  needs  to  be  addressed  is  the  specifications  of  the 
“damaged”  beam.  The  installed  beam  damage  must  produce  a  discernable 
difference  in  peak  frequencies  in  as  many  modes  as  possible  in  order  to  provide 
Equation  (37)  with  a  usable  eigenvalue  difference  vector.  In  order  to  find  a 
suitable  “damage”  condition,  numerous  FEM  models  representing  various 
candidate  conditions  must  be  developed. 

A.  MESH  PATTERNS 

Although  the  basic  geometry  of  the  beam  FEM  remains  constant,  several 
different  meshing  patterns  may  be  utilized.  The  base  requirement  is  that  a  single 
element  account  for  the  physical  alteration  made  to  represent  damage.  This  way, 
the  impact  on  the  system’s  stiffness  is  accurately  accounted  for.  Flowever,  this 
leaves  the  rest  of  the  beam  open  to  various  options  in  element  size,  number,  and 
arrangement.  To  assess  the  viability  of  these  various  mesh  arrangements,  the 
beam’s  theoretical  natural  frequencies  are  calculated.  W.D.  Pilkey’s  Formulas  for 
Stress,  Strain,  and  Structural  Matrices  provides  a  table  of  various  beam 
conditions  and  their  associated  formulae  for  calculating  natural  frequencies  for  all 
modes  [13]. 

In  order  to  replicate  these  values  as  precisely  as  possible,  a  high-fidelity 
FEM  was  generated  using  360  elements,  each  covering  0.508  cm  (0.2  in)  of  the 
beam’s  length.  As  can  be  seen  in  Table  3,  the  high-fidelity  FEM  natural 
frequencies  have  minute  differences  when  compared  to  their  theoretical 
counterparts.  Therefore,  this  mesh  pattern  is  capable  of  producing  an  accurate 
FEM.  However,  this  high-fidelity  arrangement  is  limited  with  regard  to 
implementing  various  damage  conditions.  With  a  constant  element  size,  beam 
alterations  to  represent  damage  must  ideally  be  the  same  length  as  the  element. 
An  alteration  that  is  smaller  or  larger  than  0.508  cm  (0.2  in)  would  require 
resizing  every  element  to  match  its  dimension.  Additionally,  the  high  number  of 
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elements  results  in  longer  computation  times.  While  the  element  count  (360 
elements  in  a  very  simple  FEM)  in  and  of  itself  is  not  concerning,  one  must 
understand  that  the  model  is  used  to  simulate  numerous  damaged  conditions. 
Regenerating  the  FEM  for  each  condition  certainly  has  the  potential  to  make  the 
process  time-consuming. 

As  in  alternative,  a  simple  FEM  using  the  minimum  number  of  elements 
was  also  produced.  Two  elements  were  used  to  represent  the  two  sections 
spanning  between  the  edge  of  the  beam  and  an  accelerometer  mount.  Thirteen 
additional  elements  covered  the  unaltered  beam  sections  in  between  the 
accelerometer  locations.  A  sixteenth  element  of  variable  size  was  positioned  in 
between  two  accelerometer  locations  to  cover  the  alteration  (“damage”).  The 
remaining  seventeenth  and  eighteenth  element  covered  the  remaining  space 
between  the  accelerometers  not  affected  by  the  alteration.  Figure  15  shows  the 
arrangement  of  nodes  with  the  accelerometer  locations. 


Figure  15.  Simple  FEM  Nodes  (Crosses)  and  Accelerometers  (Circles) 

Obviously,  the  minimum  number  of  elements  is  advantageous  with  regard  to 
computation  time.  Flowever,  as  seen  in  Table  3,  there  are  detectable  differences 
between  the  natural  frequencies  calculated  from  the  simple  FEM  and  the 
corresponding  theoretical  values. 

A  third  option  was  produced  using  a  transition  mesh.  This  FEM  is  similar 
to  the  simple  model,  in  that  the  beam  alteration  is  covered  by  a  variable  sized 
element.  Flowever,  in  this  case,  the  elements  surrounding  the  damage  element 
gradually  increase  in  size,  while  still  maintaining  nodes  at  accelerometer 
locations.  Figures  1 6  and  1 7  show  the  overall  mesh  and  the  transition  in  detail. 
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Figure  16.  Transition  FEM  Nodes  (Crosses)  and  Accelerometers 

(Circles) 
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Figure  17.  Detail  of  Transition  Mesh  Nodes  (Crosses) 

This  model  offers  flexibility  in  alteration  sizing  and  a  compromise  between 
element  count  and  computation  time.  As  per  Table  3,  the  difference  between  the 
transition  mesh  FEM  natural  frequency  and  the  theoretical  values  is  superior  to 
that  of  the  simple  model  (Appendix  C  provides  a  more  detailed  table  regarding 
FEM  comparison).  Due  to  these  advantages,  the  transition  mesh  FEM  was 
utilized  in  assessing  damage  condition  candidates. 


Table  3.  FEM  Natural  Frequency  Error 
(Compared  to  Theoretical  Values) 


Mode 

High  Fidelity  FEM 
(%  absolute  error) 

Simple  FEM 
(%  absolute  error) 

Transition  FEM 
(%  absolute  error) 

1 

0.0001 

0.0008 

0.0007 

2 

0.0000 

0.0054 

0.0044 

3 

0.0000 

0.0201 

0.0173 

4 

0.0000 

0.0597 

0.0489 

5 

0.0000 

0.1283 

0.1005 

6 

0.0000 

0.2178 

0.2076 

7 

0.0000 

0.3821 

0.3504 

8 

0.0000 

0.6817 

0.5064 

9 

0.0000 

0.9993 

0.7962 

10 

0.0000 

1.2631 

1 .0278 

11 

0.0000 

1 .7552 

1 .3744 

12 

0.0000 

2.4256 

2.2570 

1  through  12 

0.0001 

7.9392 

6.6918 
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B.  DAMAGE  CONDITION  SELECTION 


With  a  flexible  FEM  with  reasonable  accuracy  and  computation  time 
available,  various  damage  conditions  can  be  assessed.  Since  the  primary 
objective  of  this  research  is  experimental  validation,  a  relatively  simple  scenario 
featuring  a  single  damaged  location  is  utilized. 

The  preferred  beam  alteration  to  represent  damage  is  two  milled  notches 
across  the  width  of  the  beam.  The  simple  geometry  of  these  notches  permits 
accurate  simulation  by  the  FEM  damage  element.  The  reduced  thickness  of  the 
beam  covered  by  this  element  reduces  the  second  moment  of  inertia,  and, 
therefore,  lowers  the  element  stiffness.  Having  a  notch  on  either  side  of  the  beam 
permits  similar  beam  flexure  in  either  direction,  thus  eliminating  non-linearities  in 
the  dynamic  response.  Figures  18  and  19  show  the  milled  notches  used  to 
represent  damage  in  the  test  beam.  It  is  also  important  to  note  that  the  notch 
geometry  must  preclude  the  possibility  of  the  notch  sides  coming  in  contact  with 
each  other  as  the  beam  bends.  Such  an  occurrence  could  introduce  nonlinear 
gap  closing  in  the  beam’s  dynamic  response. 
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Figure  19.  Isometric  View  of  Milled  Notches  (“Damage”) 


With  the  acceptance  of  basic  geometry  of  the  “damage,”  various  locations 
and  notch  dimensions  have  to  be  assessed.  As  previously  stated,  the  objective  is 
to  select  a  damage  scenario  that  results  in  discernable  natural  frequency 
changes  across  as  many  modes  as  possible  (Figures  20  and  21  show  the  effect 
of  stiffness  reduction  on  natural  frequencies).  Such  a  condition  would  provide  the 
most  useful  data  via  the  eigenvalue  difference  vector  in  Equation  (37). 
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Stiffness  Reduction  at  Damaged  Element  (%) 


Figure  20.  Natural  Frequencies  as  Stiffness  Is  Reduced  at 

Damage  Location 


Figure  21 .  Change  in  Natural  Frequencies  as  Stiffness  Is  Reduced 
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A  variety  of  notch  widths,  depths,  and  locations  were  simulated. 
Specifically,  twelve  potential  elements,  five  notch  widths,  and  five  notch  depths 
were  evaluated  for  a  total  of  300  different  damage  scenarios.  Each  resultant 
change  in  natural  frequency  within  the  first  twelve  modes  was  graded  based  on 
their  detectability.  An  ideal  frequency  delta  is  detectable  by  instrumentation  with 
a  frequency  resolution  of  0.625  Hz.  An  acceptable  delta  is  detectable  with  a 
resolution  of  0.3125  Hz.  Table  4  shows  the  resultant  shift  in  natural  frequencies 
due  to  the  installation  of  damage  of  various  dimensions  at  a  particular  element. 
Based  on  these  criteria,  notches  placed  at  interior  element  5  (between 
accelerometers  5  and  6)  with  a  width  of  0.476  cm  (0.1875  in)  and  a  depth  of 
0.238  cm  (0.09375  in).  This  alteration  creates  frequency  differences  in  ten  modes 
that  are  detectable  with  0.625  Hz  resolution. 


Table  4.  Natural  Frequency  Differences  for  Damage  at  Element  5 


Damage  Width  (cm) 

0.47625 

Damage  Depth  (cm) 

0.15875 

0.3175 

0.47625 

0.635 

0.79375 

Mode  1 

0.017929 

0.060279 

0.181296 

0.667838 

4.056329 

Mode  2 

0.06472 

0.226432 

0.686749 

2.339248 

8.509918 

Mode  3 

0.016528 

0.057506 

0.172764 

0.573156 

1.94215 

Mode  4 

0.066833 

0.240258 

0.74249 

2.561609 

9.170284 

Mode  5 

0.269053 

0.95183 

2.844721 

8.823834 

22.80683 

Mode  6 

0.069577 

0.244716 

0.721788 

2.164174 

5.379531 

Mode  7 

0.152492 

0.539591 

1.625137 

5.197318 

14.52639 

Mode  8 

0.65211 

2.277295 

6.615923 

18.82277 

39.89659 

Mode  9 

0.200309 

0.693144 

1.970191 

5.332921 

10.86344 

Mode  10 

0.237565 

0.832377 

2.44729 

7.264367 

17.1688 

Mode  1 1 

1.147385 

3.962996 

11.21054 

29.60283 

55.45561 

Mode  12 

0.359958 

1.22922 

3.38744 

8.496358 

15.57912 
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V.  ABC  SELECTION  AND  COMPOSITE  MATRIX  ASSEMBLY 


With  a  viable  damage  scenario  chosen,  potential  composite  sensitivity 
matrices  must  be  constructed  and  evaluated  before  proceeding  with  the 
experimental  portion  of  the  research.  As  previously  stated,  many  ABC 
configurations  (with  varying  ABC  types,  locations,  and  characteristics)  are 
applied  to  create  a  sensitivity  row  library.  Orthogonal  projection  is  then  employed 
to  evaluate  each  row  and  make  selections  that  minimize  matrix  linear 
dependence.  The  composite  sensitivity  matrix  is  then  utilized  in  a  simulated 
damage  detection  scenario  to  gauge  its  performance.  Simulation  conditions, 
featuring  damage  of  varying  magnitude  and  at  multiple  locations,  is  used  for  each 
composite  sensitivity  matrix  in  order  to  provide  a  consistent  assessment. 

A.  BASELINE  SENSITIVITY  MATRIX 

To  act  as  a  control  for  sensitivity  library  comparison,  the  baseline  (no 
ABCs)  sensitivity  matrix  of  fourteen  rows  and  fourteen  columns  provides  the  first 
set  of  sensitivity  data  to  be  assessed  via  the  orthogonal  projection  method.  The 
MDI  and  PDI  for  the  matrix  is  shown  in  Figure  22. 


MDI 


1  2  3  4  5  6  7  8  9  10  11  12  13  14 

Row  Number 


PDI 


1  2  3  4  5  6  7  8  9  10  11  12  13  14 

Column  Number 


Figure  22.  Baseline  [S]  MDI  and  PDI 
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When  subjected  to  the  damage  detection  simulation,  the  “detected”  damage 
results  are  completely  inaccurate  when  compared  to  the  “true”  damage  locations 
and  magnitudes  (see  Figure  23).  All  quantitative  results  are  summarized  in  Table 
5. 
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Figure  23.  Simulated  Damage  Detection  Results  with  Baseline  [S] 

B.  PIN-BASED  ABC  LIBRARIES 

Pin-based  ABCs  are  simple  to  apply  from  a  coding  perspective  and  are 
relatively  quick  to  assess,  as  they  do  not  have  adjustable  characteristics  (unlike  a 
spring,  which  can  have  varying  stiffness).  The  first  ABC-augmented  library  to  be 
evaluated  consists  of  matrices  produced  from  the  baseline  configuration  (no 
ABCs  applied)  and  single  pin  scenarios.  For  each  scenario,  one  pin  is  applied  at 
a  translational  DOF  associated  with  an  accelerometer  location  on  the  test  article. 
In  total,  this  amounts  to  sixteen  ABC  configurations  (one  baseline  and  fifteen 
single  pin  configurations  for  each  accelerometer  location).  By  arranging  these 
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matrices  in  accordance  with  Equation  (35),  a  sensitivity  library  of  224  rows  and 
fourteen  columns  is  formed.  From  the  library,  orthogonal  projection  only  selected 
rows  created  by  single  pin  configurations.  The  selected  rows  from  this  expanded 
library  resulted  in  a  sharp  drop  in  composite  matrix  linear  dependence,  as 
indicated  in  Figure  24.  This,  in  turn,  provided  remarkably  improved  damage 
detection  results,  as  per  Figure  25. 


MDI 


Row  Number 


PDI 


Column  Number 


Figure  24.  Single  Pin  Composite  Sensitivity  Matrix  MDI  and  PDI 
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Figure  25.  Simulated  Damage  Results  with  Single  Pin 

Composite  Matrix 

The  current  sensitivity  library  featuring  the  baseline  matrix  and  all  single 
pin  scenarios  can  be  further  expanded  by  adding  rows  generated  from  two  pin 
configurations.  The  number  of  available  sensitivity  rows  grows  extensively  to  a 
count  of  1694.  It  is  noteworthy  that  orthogonal  projection  only  selected  rows  from 
two  pin  ABC  configurations.  The  greater  number  of  candidate  rows  allows  for  a 
further  decrease  the  matrix’s  linear  dependence  even  further,  as  evidenced  by 
the  associated  MDI  and  PDI  (see  Figure  26).  Once  again,  simulated  damage 
detection  results  are  promising  (see  Figure  27). 
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Figure  26.  Two  Pin  Composite  Sensitivity  Matrix  MDI  and  PDI 


Element  No. 


Figure  27.  Simulated  Damage  Results  with  Two  Pin  Composite  Matrix 

Adding  sensitivity  rows  generated  from  three  pin  ABC  configurations 
exponentially  grows  the  size  of  the  sensitivity  library  to  8064  rows.  In  a  manner, 
similar  to  the  previous  cases,  orthogonal  projected  selected  rows  from 
configurations  with  the  highest  pin  count.  Linear  dependency  is  further  reduced, 
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but  the  improvement  is  not  as  significant  compared  to  previous  additions.  MDI 
and  PDI  barely  show  the  reduction  in  linear  dependency  (Figure  28)  and  the 
damage  simulation  exhibits  minutely  improved  results  (Figure  29).  All  quantitative 
results  associated  with  pin-based  ABCs  are  found  in  Table  5. 


PDI 


1  2  3  4  5  6  7  8  9  10  11  12  13  14 

Column  Number 


Figure  28.  Three  Pin  Composite  S  Matrix  MDI  and  PDI 


Figure  29.  Simulated  Damage  Results  with  Three  Pin  Composite  Matrix 
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C.  ALTERNATIVE  ABC  TYPES 

As  previously  alluded,  there  are  additional  ABCs  besides  pins.  R.J.C. 
Konze  conducted  research  on  the  use  of  springs  and  point  masses  as  ABCs. 
Based  on  conclusion  and  recommendations  from  “Synthetic  modification  in  the 
frequency  domain  for  finite  element  model  update  and  damage  detection,”  a 
single  spring  with  stiffness  ranging  from  175.13  N/mm  (1000  Ib/in)  to  1751.3  N/ 
mm  (10000  Ib/in)  generates  the  most  useful  sensitivity  rows  for  a  free-free  beam. 
For  point  masses,  the  recommended  values  range  from  0.45359  kg  (1  Ibm)  to 
1.8143  kg  (4  Ibm)  [12]. 

A  separate  sensitivity  library  was  created  with  row  inputs  from  the  baseline 
sensitivity  matrix,  single  pin  ABCs,  single  spring  ABCs,  and  single  point  mass 
ABCs.  As  with  the  previous  cases,  each  ABC  was  applied  to  each  node 
corresponding  to  the  fifteen  accelerometer  locations.  The  library  size  reached 
3164  rows.  Interestingly,  of  the  fourteen  composite  sensitivity  matrix  rows, 
orthogonal  projection  selected  twelve  rows  from  spring-based  ABCs,  two  from 
pin-based  ABCs,  and  none  from  point  mass-based  ABCs.  Compared  to  the 
composite  matrix  constructed  from  a  single  pin  ABCs,  alone,  the  spring-pin 
composite  exhibited  lower  linear  dependency  (see  Figure  30).  This  composite 
sensitivity  matrix  also  produced  satisfactory  simulated  damage  detection  results 
(see  Figure  31).  Quantitative  results  are  summarized  in  Table  5. 
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Figure  30.  Single  Pin,  Spring,  Point  Mass  Composite  Matrix 

MDI  and  PDI 


Figure  31 .  Simulated  Damage  Results  with  Single  Pin,  Spring, 

Point  Mass  Composite  Matrix 

D.  CONCLUDING  REMARKS 

In  total  five  different  composite  sensitivity  matrices  were  constructed  and 

evaluated  (including  the  baseline  sensitivity  matrix).  The  specific  rows  used  in 

their  construction  are  listed  in  Appendix  C.  A  quantitative  comparison  is  shown  in 
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Table  5.  Maximum  false  damage  error  is  the  value  of  the  largest  damage 
detection  return  for  a  non-damaged  element.  Maximum  damaged  element  error 
is  the  largest  difference  between  the  “true”  damage  value  and  the  “detected” 
value  for  a  damaged  element.  Total  stiffness  error  is  simply  the  sum  of  all 
differences  between  “true”  and  “detected”  damage  across  all  elements. 


Table  5.  Composite  Sensitivity  Matrix  Comparison 


Available 

ABCs 

Library 

Size  (rows) 

Condition 

Number 

Max  false 
damage 
error  (%) 

Max 

damaged 
element 
error  (%) 

Total 
stiffness 
error  (%) 

None 

14 

13786 

119.00 

113.11 

852.49 

1  Pin 

224 

11.204 

0.17028 

0.50788 

1 .6279 

1  Pin, 

2  Pins 

1694 

6.9989 

0.14012 

0.33576 

1.2594 

1  Pin, 

2  Pins, 

3  Pins 

8064 

5.7334 

0.14128 

0.22850 

1 .0642 

1  Pin, 

1  Spring, 

1  Pt.  Mass 

3164 

9.3595 

0.43835 

0.45066 

2.5892 

It  is  obvious  that  the  use  of  ABCs  to  improve  matrix  conditioning  is 
required.  It  is  also  clear  that  increasing  the  number  of  pins  further  improves 
condition  number  and  reduces  overall  error.  However,  there  is  a  diminishing 
return  on  increasing  the  pin  count.  As  more  pins  are  added,  the  accuracy 
improves  by  smaller  magnitudes.  On  another  note,  the  implementation  of  spring- 
ABCs  was  deemed  worthwhile,  per  orthogonal  projection  selection.  Compared  to 
a  single  pin  composite  matrix,  the  pin-spring  composite  matrix  featured  a  lower 
condition  number.  While,  the  total  stiffness  error  was  slightly  higher  for  the  pin¬ 
spring  composite,  it  did  have  a  lower  maximum  damaged  element  error. 
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VI.  FEM  UPDATING 


With  several  viable  composite  sensitivity  matrices  available,  all 
prerequisites  for  model  updating  are  met.  The  FEM  to  be  updated  is  the  same  as 
the  Simple  FEM  described  in  Chapter  IV,  minus  the  dedicated  element  for 
damage  insertion.  This  model  provides  fourteen  interior  elements  representing 
the  sections  of  the  beam  encompassed  by  the  mounted  accelerometers.  These 
are  the  elements  that  will  be  updated  during  the  model  update  process.  The 
remaining  part  of  the  beam  (the  space  between  the  beam’s  edges  and  the 
outboard  accelerometers)  are  represented  by  two  additional  elements  that  will 
remain  constant. 


In  an  effort  to  bring  FEM  natural  frequencies  into  close  alignment  with  the 
experimentally  obtained  values,  the  following  multiple-iteration  update  process  is 
used: 


1. 


Let{M  =  {V'’}  [?]=pw 


,  {AA}=|AA(,)|  (first  iteration  /=1) 


2.  Update  the  FEM  with 


A p 


(0 


3. 


Recalculate  the  sensitivity  rows 
calculated  in  Step  2. 


and  natural  frequencies 


4.  Construct  ‘  >  using  the  newly  calculated  natural  frequencies 

along  with  the  frequencies  already  found  from  curve-fitting  the 
synthesized  FRF. 


5. 


Solve  for 


6. 


If 


{a/+1)}  Ja/)} 


,  End. 
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7. 


A 

Otherwise,  let 
and  return  to  Step  1 . 


Four  separate  modal  tests  were  conducted  on  the  intact  test  article,  using 
a  variety  of  FRF  frequency  resolutions.  The  first  test,  using  the  single  impact 
method,  features  a  resolution  of  0.3125  Hz.  The  second  test,  also  using  single 
impact,  has  a  resolution  of  0.08  Hz.  However,  the  requisite  response  recording 
time  for  this  resolution  is  twelve  seconds.  As  this  relatively  long  recording  time 
makes  this  method  increasingly  susceptible  to  noise,  an  additional  test  at  0.08  Hz 
resolution  was  conducted  using  the  random  impact  method.  All  utilized  FRFs  at 
this  resolution  were  from  the  random  impact  test.  A  fourth  test  at  a  resolution  of 
0.02  Hz  was  performed  with  random  impact.  The  measured  FRFs  are  shown  in 
Figure  32. 

It  is  important  to  note  that  the  peak  frequencies  are  affected  by  damping 
and  are,  therefore,  not  the  natural  frequencies.  Due  to  the  free-free  condition  of 
the  test  article,  damping  is  kept  to  a  minimum.  In  order  to  obtain  the  most 
accurate  results,  a  single-DOF  curve  fitter  is  applied  to  peaks  of  interest.  This 
algorithm,  developed  A.  M.  Rinawi  and  R.  W.  Clough,  analyzes  the  peak, 
estimates  the  damping  ratio,  and  provides  the  natural  frequency  [14].  These 
values  can  then  be  compared  to  the  associated  FEM  values.  The  differences 
between  these  frequencies  can  then  be  used  to  solve  for  the  desired  flexural 
rigidity  adjustments  in  accordance  with  Equation  (37). 
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Figure  32.  Experimentally  Obtained  FRFs  from 

Each  Modal  Test 

A.  COMPLICATIONS  REGARDING  MODEL  UPDATE 

In  an  effort  to  obtain  the  best  experimental  results,  a  modal  test  using  the 
random  impact  method  with  a  high  frequency  resolution  of  0.02  Hz  was 
conducted.  The  experimentally  obtained  FRF  is  shown  in  Figure  33.  With  the  H 
matrix  obtained,  the  ABCs  can  be  applied  to  the  experimental  data  using 
synthesis.  The  first  ABC  sets  to  be  applied  are  those  used  to  construct  the  two- 
pin  composite  sensitivity  matrix.  As  previously  discussed  in  Chapter  V,  this  matrix 
featured  a  very  low  condition  number  while  still  maintaining  a  reasonably  sized 
sensitivity  library  (keeping  computation  times  reasonable).  Based  on  the  merits  of 
high  FRF  frequency  resolution  and  ABCs  associated  with  minimal  sensitivity 
linear  dependence,  an  exceptional  model  updating  result  is  expected. 

However,  when  artificial  pins  were  applied  at  accelerometer  locations  1 
and  7,  the  resultant  H*  matrix  produced  the  accompanying  FRF  shown  in  Figure 
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33.  This  synthesized  FRF  is  severely  distorted  compared  to  its  baseline 
counterpart. 


Figure  33.  FRF  Distortion  with  0.02  Hz  Resolution  and 

Two  Pin  ABCs 

Other  ABC  configurations  required  by  the  composite  sensitivity  matrix  also 
produce  varying  degrees  of  distortion.  Figure  34  shows  the  synthesized  FRF 
from  a  different  ABC  set.  This  figure  highlights  the  influence  of  artificial  pin 
location  on  FRF  synthesis. 
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Figure  34.  FRF  Distortion  with  0.02  Hz  Resolution  and 

Two  Pin  ABCs 

While  Figure  34  shows  a  less  severe  distortion  case,  in  general,  the  extent  of 
distortion  makes  peak  identification  impossible.  This  precludes  the  identification 
of  their  respective  frequencies.  Without  the  peak  frequencies,  the  experimental 
natural  frequencies  remain  unattainable,  flexural  rigidity  adjustments  remain 
unsolved,  and  FEM  updating  becomes  an  impossible  task. 

Additional  attempts  at  FRF  synthesis  are  made  using  lower  resolutions  in 
an  effort  to  avoid  this  debilitating  distortion.  Figures  35  through  37  show  samples 
of  FRFs  synthesized  from  a  variety  of  ABCs  and  resolutions. 
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Figure  35. 
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Figure  37.  Synthesized  FRF  with  Two  Pins  at 
0.08  Hz  Resolution 

Based  on  these  FRFs,  it  is  clear  that  distortion  is  a  function  of  both  ABC  count 
and  resolution.  Figures  35  and  36  show  increased  distortion  due  to  using  two 
pins  instead  of  one.  Figures  36  and  37  show  increased  distortion  due  to  higher 
resolution. 

With  viable  synthesized  FRFs  obtained,  model  update  is  first  attempted 
using  single  pin  ABCs  and  a  resolution  of  0.3125  Hz.  Due  to  its  minimal 
distortion,  peaks  and  their  frequencies  are  readily  identified.  Due  to  the  minimal 
damping  of  the  test  article,  these  peak  frequencies  are  extremely  close  to  the 
actual  natural  frequencies.  With  the  experimentally  obtained  natural  frequency 
inputs,  the  flexural  rigidity  adjustments  are  computed.  Each  adjustment  can  then 
be  individually  applied  to  its  respective  interior  element. 

It  is  important  to  note  that  E.  Damanakis’s  research  on  model  updating 
advocated  averaging  all  flexural  rigidity  adjustments  and  uniformly  applying  the 
average  to  the  FEM  [15].  While  this  research  presents  a  successful  case,  there  is 

a  critical  caveat.  This  approach  requires  average  values  that  are  positive  be 
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completely  disregarded.  Uniformly  applying  a  positive  flexural  rigidity  adjustment 
would  result  in  higher  FEM  natural  frequencies,  thus  producing  larger  eigenvalue 
differences  between  the  model  and  the  test  article  [15].  The  fact  that  positive 
flexural  rigidity  adjustment  averages  prevent  any  form  of  model  update  makes 
this  approach  far  from  ideal.  Figure  38  highlights  this  deficiency.  It  is  for  this 
reason  that  individual  element  adjustments  are  utilized  here. 


Figure  38.  Example  of  Disregarded  Flexural  Rigidity 
Adjustments.  Source:  [15].  Their  positive  average  (green  line) 
precluded  their  consideration. 


As  previously  discussed,  this  model  updating  process  employs  a  multiple 
iteration  procedure.  With  every  iteration,  as  the  FEM  frequencies  converge  on 
the  experimental  values,  the  required  adjustments  should  decrease.  Figure  39 
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exhibits  the  calculated  adjustments  over  the  course  of  two  model  update 
iterations. 


Figure  39.  Model  Update  with  Single  Pin  ABC  at 

0.3125  Hz  Resolution.  All  modes  were  available  from 
sensitivity  library. 

This  figure  clearly  shows  the  flexural  rigidity  adjustments  growing  in  magnitude, 
indicating  divergence.  This  unfortunate  result  is  another  example  of  a  failed 
model  update. 

The  cause  of  this  phenomenon  is  found  in  the  lower  modes  of  the 
synthesized  FRF.  The  percent  eigenvalue  differences  between  the  FEM  and  the 
experimental  data  are  noticeably  larger  in  several  instances  where  a  low  mode 
frequency  is  used.  Appendix  C  features  tables  containing  all  the  initial  eigenvalue 
differences.  The  consequence  of  these  larger  differences  is  that  they  exceed 
eigenvalue  sensitivity  linearity,  thus  resulting  in  a  failure  to  converge.  These 
abnormally  large  percent  differences  can  be  attributed  to  one  of  three  separate 
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causes,  depending  on  the  scenario.  The  first  source  of  error  is  increased  peak 
distortion  at  the  lower  modes  of  the  synthesized  FRF.  This  trend  is  visible  in 
Figure  37  and  an  extreme  case  shown  in  detail  in  Figure  40. 


A/tMaM  Pn  Appfeed  at  U 


I  o. 

i 

s 


Figure  40.  Distortion  of  an  FRF  Peak  at  Mode  2. Based  on  FEM 
predictions,  the  peak  should  be  in  the  vicinity  of  19.572  Hz. 


Due  to  this  distortion,  peak  selection  is  skewed,  resulting  in  large  eigenvalue 
differences.  Another  cause  of  error  is  the  shifting  of  modes  into  extremely  low 
frequencies.  At  frequencies  of  six  Hz  or  less,  the  impact  testing  can  fail  to 
produce  an  FRF  peak  in  its  entirety.  An  example  is  shown  in  Figure  41 . 
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Figure  41 .  Example  of  Shifting  Mode  Frequencies. 

The  first  mode  shifts  to  such  a  low  frequency  that  a  peak  is  not 
generated.  A  peak  should  exist  in  the  vicinity  of  0.6445  Hz  (per 

initial  FEM). 


A  third  source  of  error  is  more  subtle.  There  are  instances  where  a  peak  is 
produced  with  minimal  distortion,  yet  still  produces  large  eigenvalue  differences. 
One  such  case  is  shown  in  Figure  42. 
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Artificial  Pin  Applied  at  Location  9 
No  Artificial  Pins  Applied | 


Frequency (Hz) 


Figure  42.  Clean  Peak  at  Mode  1  at  10.625  Hz.  Associated  FEM 

frequency  is  9.552  Hz.  Resultant  eigenvalue  percent  difference 

is  22.362%. 


Such  large  percent  eigenvalue  differences,  like  the  example  in  Figure  42,  stand 
out  from  the  rest  of  eigenvector.  All  of  the  initial  percent  differences  are 
tabularized  in  Table  11  in  Appendix  C.  Referring  to  Table  10  in  conjunction  with 
Table  11,  it  is  clear  that  the  larger  differences  are  found  in  the  lowest  modes. 
This  eigenvalue  difference  can  be  further  exacerbated  by  the  frequency 
resolution,  which  has  a  larger  effect  on  low  frequency  modes  as  opposed  to  high 
modes.  Differences  of  such  magnitude  can  easily  exceed  sensitivity  linearity. 

It  is  clear  that  low  modes  introduce  potential  error  sources  in  both  modal 
testing  and  ABC  synthesis.  These  errors  can  lead  to  high  percent  eigenvalue 
differences  exceeding  sensitivity  linearity,  and,  in  turn,  divergence  during  model 
update.  However,  the  lowest  modes  are  not  the  only  potential  source  of 
divergence.  As  previously  shown  in  Table  3,  FEM  generated  natural  frequencies 
begin  to  deviate  from  theoretical  values,  especially  when  the  FEM  has  a 
relatively  low  element  count.  This  is  highlighted  in  Figure  43. 
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Figure  43.  Deviation  of  FEM  Natural  Frequencies  from  Theoretical 

Values  at  Higher  Modes 


Fortunately,  none  of  the  composite  sensitivity  matrices  used  in  model  updating 
utilize  modes  higher  than  10.  This  mode  does  not  feature  significant  natural 
frequency  deviation.  However,  this  phenomenon  should  be  taken  into 
consideration. 

B.  MODEL  UPDATE  RESULTS 

Based  on  these  potential  error  sources,  a  conservative  approach  is  taken, 
and  the  sensitivity  library  is  limited  to  rows  utilizing  modes  4  through  8  only  for 
future  model  update  attempts.  Since  the  natural  frequency  deviation  at  higher 
modes  is  not  extreme,  as  shown  by  Figure  43,  an  additional  model  update 
scenario  using  a  library  with  modes  4-14  is  also  tried  for  the  sake  of  comparison. 
Scenarios  where  multiple  pins  are  synthesized  into  high  resolution  FRFs  are 
avoided  entirely. 
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These  new  model  update  attempts  result  in  remarkable  success.  These 
successful  model  updating  cases  and,  and  their  conditions,  are  summarized  in 
Table  6. 


Table  6.  Updated  FEM  Case  Comparison 


Updated  FEM 
Case 

Resolution 

(Hz) 

ABCs 

Available 

Modes 

Test  Method 

A 

0.3125 

1  Pin 

4-14 

Single  Impact 

B 

0.3125 

1  Pin 

4-8 

Single  Impact 

C 

0.3125 

2  Pins 

4-8 

Single  Impact 

D 

0.08 

1  Pin 

4-8 

Random 

Impact 

E 

0.08 

2  Pins 

4-8 

Random 

Impact 

F 

0.02 

1  Pin 

4-8 

Random 

Impact 

The  first  encouraging  indication  is  the  reduction  in  flexural  rigidity 
adjustments  with  each  successive  update  iteration.  Figures  44  through  49  exhibit 
this  trend. 
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Figure  44.  FEM  Case  A  Flexural  Rigidity  Updates  Over  Three  Iterations 
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Figure  45.  FEM  Case  B  Flexural  Rigidity  Updates  Over  Three  Iterations 
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Figure  46.  FEM  Case  C  Flexural  Rigidity  Updates  Over  Three  Iterations 
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Figure  47.  FEM  Case  D  Flexural  Rigidity  Updates  Over  Three  Iterations 
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Stiffness  Adjustment  (Percent  of  Nominal  Stiffness) 


Figure  48.  FEM  Case  E  Flexural  Rigidity  Updates  Over  Three  Iterations 
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Figure  49.  FEM  Case  F  Flexural  Rigidity  Updates  Over  Three  Iterations 
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With  these  updates  applied  individually  to  their  respective  interior  elements,  the 
percent  differences  between  FEM  natural  frequencies  and  their  experimental 
counterparts  are  significantly  reduced.  Figure  50  shows  the  extent  of  the  model 
updating  success.  This  plot  features  the  initial  FEM,  all  six  updated  FEM  cases, 
and  a  manually  adjusted  FEM  implemented  by  Damanakis  in  his  research 
(“Adjusted  FEM”  in  Figure  50). 


-*■- Test-Updated  FEM  C  Test-Updated  FEM  D  — t— ' Test-Updated  FEM  E  - Test-Updated  FEM  F 

Figure  50.  Reduction  of  Eigenvalue  Percent  Differences  between  FEMs 

and  Experimental  Values 

C.  CONCLUDING  REMARKS 

As  seen  is  this  chapter,  combining  multiple  pin  ABCs  with  high  resolution 
FRFs  produce  severe  distortion,  and  such  combinations  must  be  avoided.  The 
lowest  modes  introduce  multiple  error  sources  that  prevent  converging  model 
updates,  and  their  use  should  also  be  avoided.  While  FEM  natural  frequencies 
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deviate  from  theoretical  values  at  higher  modes,  these  differences  are  not  large 
enough  to  preclude  successful  model  updating. 

Per  Figure  50,  a  three-iteration,  individual  element  model  updating 
approach  is  superior  to  a  uniformly  applied,  manual  FEM  adjustment.  It  is 
important  to  note  that  updating  cases  utilizing  two  artificial  pins  resulted  in  less 
difference  reduction  than  single  pin  cases.  These  relatively  mediocre  results 
might  offset  any  advantage  the  two  pins  provide  in  terms  of  sensitivity  condition 
number. 
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VII.  DAMAGE  DETECTION 


With  multiple  updated  FEMs  available,  damage  detection  testing  can  be 
conducted.  As  previously  discussed,  this  involves  comparing  the  natural 
frequencies  of  the  updated  models  to  the  experimentally  obtained  natural 
frequencies  of  the  “damaged”  test  article.  Of  note,  the  damage  detection  case 
notation  corresponds  to  the  model  updating  cases. 

A.  DAMAGE  DETECTION  COMPLICATIONS 

The  first  attempts  at  damage  detection  utilize  updated  FEMs  C  and  E, 
which  feature  two  pin  ABCs.  As  mentioned  in  Chapter  VI,  two  pin  model  updating 
was  not  as  effective  as  the  cases  using  single  pins.  Flowever,  as  shown  in 
Chapter  V,  two  pin  ABCs  provided  exceptional  damage  detection  results  in 
simulation.  The  resultant  damage  detection  graphs  for  cases  C  and  E  are  shown 
in  Figures  51  and  52.  The  damage  should  be  manifested  as  a  distinct  drop  in 
flexural  rigidity  at  the  correct  damage  location  element  (interior  Element  5). 
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Figure  51 .  Case  C  Damage  Detection  Results.  The  correct  damage 
location  is  shown  at  Element  5,  but  with  significant  false  alarms. 


Element  No. 


Figure  52.  Case  E  Damage  Detection  Results.  The  correct  damage 
location  at  is  shown  at  Element  5,  but  with  significant  false 

alarms. 
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Both  Cases  C  and  E  correctly  locate  the  test  article’s  installed  “damage,”  as 
shown  with  the  large  reduction  in  flexural  rigidity  at  Element  5.  However,  there 
are  numerous  false  alarms  with  significant  magnitudes.  These  false  alarms  are 
large  enough  to  completely  compromise  the  ability  to  distinguish  damaged 
elements  from  undamaged  ones.  Table  7  summarizes  the  quantitative  damage 
detection  results. 

An  analysis  of  these  disappointing  results  quickly  reveals  the  limiting 
factor.  As  with  model  updating,  even  mild  FRF  peak  distortion  can  have  a 
negative  impact  on  damage  detection,  even  when  modes  1  through  3  are 
excluded.  Figure  53  shows  the  effect  of  two  pin  ABC  configurations  on  the 
synthesized  FRF. 


Artificial  Pin  Applied  at  Locations  12  and  13 
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Figure  53.  Detail  of  FRF  Peak  Distortion  with  Two  Pin  ABC 

Configuration 
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B.  DAMAGE  DETECTION  RESULTS 


Due  to  the  negative  impact  of  two  pin  ABC  configurations  on  the 
synthesized  FRF,  further  damage  detection  efforts  are  limited  to  single  pin 
configurations.  Damage  detection  results  from  Cases  A,  B,  D,  and  E  are 
exhibited  in  Figures  54  through  57. 
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Figure  54.  Case  A  Damage  Detection  Results  with  Clear  Damage 

Location 
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Change  in  Stiffness  (Percent  of  Nominal  Stiffness) 
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Figure  55.  Case  B  Damage  Detection  Results  with  Clear  Damage 

Location 


Figure  56.  Case  D  Damage  Detection  Results  with  Clear  Damage 

Location 
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Figure  57.  Case  F  Damage  Detection  Results  with  Clear  Damage 

Location 


Based  on  these  figures,  it  is  clear  that  the  damaged  element  can  be  readily  and 
consistently  identified.  The  quantitative  results  are  summarized  in  Table  7. 

C.  CONCLUDING  REMARKS 

With  regard  to  damage  detection,  the  advantages  of  single  pin  ABC 
configurations  over  two  pin  configurations  is  obvious.  Table  7  compares  the 
damage  detection  performance  of  each  of  the  test  cases.  Total  noise  indicates 
the  sum  of  all  returns  associated  with  an  undamaged  element.  Maximum  false 
detection  is  the  largest  single  false  damage  return  on  an  intact  element. 
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Table  7.  Damage  Detection  Performance  Comparison 


Damage  Detection 
Case 

Total  Noise  (%) 

Maximum  False 
Detection  (%) 

A 

44.879 

8.3842 

B 

32.425 

4.9576 

C 

114.14 

19.489 

D 

28.953 

5.9607 

E 

97.431 

17.234 

F 

29.357 

4.2635 

Table  7  emphasizes  the  superior  results  produced  by  single  pin  ABC 
configurations.  It  is  also  important  to  note  that  Case  A  has  slightly  higher  error 
than  the  other  single  pin  cases.  However,  the  calculated  flexural  rigidity  drop  at 
the  damage  element  has  a  higher  magnitude  than  the  other  cases,  and, 
therefore,  maintains  a  similar  prominence  against  the  false  returns. 
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VIII.  CONCLUSION 


A.  KEY  OUTCOMES 

The  first  key  area  of  research  involved  the  assessment  of  various  FEM 
meshing  schemes  to  simulate  potential  damage  scenarios.  The  concept  of  a 
transitional  mesh  in  which  mesh  size  gradually  expands  with  increased  distance 
from  installed  damage  offered  a  compromise  of  natural  frequency  accuracy  and 
reasonable  computation  time.  Additionally,  it  provides  the  flexibility  required  to 
assess  damage  of  various  dimensions.  While  the  savings  in  computation  time  is 
not  critical  for  a  simple  beam,  it  will  become  increasingly  important  as  test  articles 
of  increasingly  complex  geometry  are  tested. 

With  regard  to  using  ABCs  to  reduce  linear  dependency  in  the  composite 
sensitivity  matrix,  it  is  obvious  that  there  are  many  various  combinations  and 
configurations  that  can  be  employed.  Using  orthogonal  projection  to  select  the 
sensitivity  rows,  it  is  clear  that  multiple  ABCs  are  ideal  for  driving  condition 
number  toward  unity.  When  comparing  single  pin  to  double  pin  to  triple  pin 
configurations,  orthogonal  projection  selects  the  rows  generated  from  triple  pin 
conditions.  As  shown  with  damage  detection  simulations,  the  sensitivity  matrices 
with  lower  condition  numbers  have  superior  performance  with  regard  to  detection 
accuracy.  However,  using  multi-pin  combinations  greatly  increases  the  size  of 
the  sensitivity  library  and  results  in  longer  computation  times.  With  every 
additional  pin  added,  there  is  a  diminishing  return  on  performance  improvement. 

The  comparison  of  pin-based  ABCs  to  other  types,  such  as  spring  and 
point  masses,  reveals  interesting  results.  When  provided  with  a  library  with  single 
pin,  single  spring,  and  single  point  mass  ABCs,  orthogonal  projection  constructed 
a  composite  matrix  made  from  mostly  spring-based  rows  and  several-pin  based 
rows.  While  overall  accuracy  was  not  improved  in  comparison  to  the  single-pin 
matrix,  the  mixed  ABC  matrix  offered  reduced  maximum  error  in  the  damaged 
elements. 
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Turning  to  the  experimental  portion  of  the  research,  it  quickly  became 
clear  that  raw  FRF  data  is  susceptible  to  distortion.  The  magnitude  of  distortion  is 
a  function  of  FRF  frequency  resolution  and  the  number  of  ABCs  synthesized  into 
the  /-/matrix.  Higher  ABC  counts  and  higher  resolution  lead  to  higher  distortion. 
Lower  modes  are  particularly  vulnerable  to  this  phenomenon,  as  well  as  sources 
of  experimental  error. 

These  complications  can  be  avoided  through  the  judicious  use  of  select 
ABCs,  FRF  resolution,  and  modes.  With  the  proper  conditions  employed, 
iterative  model  updating  can  be  employed.  This  allows  for  flexural  rigidity 
adjustments  to  be  applied  to  every  individual  interior  element.  When  compared  to 
manual  adjustments  uniformly  applied  to  the  entire  FEM  (as  performed  in 
Damanakis’s  research),  the  iterative  model  update  is  shown  to  be  superior. 

The  concerns  regarding  distortion  and  experimental  error  in  model 
updating  also  apply  for  damage  detection.  Even  mild  FRF  peak  distortion  can 
compromise  damage  detection.  Even  if  damage  is  correctly  located  via 
indications  of  flexural  rigidity  reduction,  false  returns  can  compromise  the 
process’s  effectiveness.  However,  as  with  model  updating,  prudent  usage  of 
ABCs,  resolution,  and  modes  can  lead  to  clear,  consistent  damage  detection. 

B.  RECOMMENDATIONS 

Based  on  obtaining  these  promising  results,  the  next  stage  of  research 
should  be  optimization.  Additional  ABC  combinations  should  be  investigated  in 
an  effort  to  drive  sensitivity  matrix  condition  number  toward  unity.  Since  this  will 
greatly  expand  the  sensitivity  library,  high-power  computers  should  be  employed 
to  avoid  unreasonable  computation  times.  On  the  experimental  side,  various 
curve  fitters  should  be  evaluated  to  generate  an  FRF  that  is  less  susceptible  to 
distortion.  Additionally,  the  existing  MATLAB  code  should  be  optimized. 
Currently,  existing  code  requires  a  large  amount  of  user  input.  Any  automation 
that  can  be  integrated  into  the  program  would  be  a  welcome  improvement.  This 
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improved  model  updating  and  damage  detection  ability  should  be  evaluated  with 
a  beam  featuring  multiple  damage  locations  of  various  dimensions. 

With  an  optimized  process,  more  realistic  scenarios  can  be  tested.  Instead 
of  milling  notches  into  the  beam,  cracking  should  be  induced  in  multiple 
locations.  Another  possible  test  would  utilize  a  heavily  corroded  test  article. 
Additional  features  can  be  evaluated  such  a  welds  and  bolted  connections.  With 
these  situations  investigated,  test  articles  with  more  complex  geometry  can  be 
used.  Instead  of  using  a  simple  rectangular  beam,  I,  T,  and  box  beams  can  be 
utilized.  Eventually,  more  sophisticated  structures  such  as  trusses  and  frames 
can  be  tested. 


Figure  58.  A  Potential  Future  Test  Article;  an  Available  Frame 
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APPENDIX  A.  TEST  ARTICLE  DETAILS 


Customer  Name 


Customer  PO# 


Shipper  No  Heat  Number 


Naval  Postgraduate  School 


KFS66964 


856510  201620673 


S, 

5  ? 


1- 


Figure  59.  Test  Beam  Inspection  Report 
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Figure  60.  Drawing  of  “Intact”  Test  Article 
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Figure  61 .  Drawing  of  “Damaged”  Test  Article 
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APPENDIX  B.  INSTRUMENTATION  DETAILS 


PCB  CALIBRATION  CERTIFICATE 

IMPULSE  FORCE  HAMMER 

086B03 


Model  No. 
Serial  No. 

Range  _ 


2195 


Customer:  \5upply 

ft>  sTtftAOu  A  tC.  SCrtiO  U 


0-500 


lb 


Linearity  error 


1L0 _ % 

2000 


Discharge  Time  Constant 
Output  Impedance  *0A  ohms 
Output  Bias  9.30  woitt 

Traceable  to  NJJS  through _ 737/236905 

Initials  f /J  Date 


Invoice  No.: 

_  35777 


.  Calibration  Specification  MIL-STO  45662 


Accelerometer: 
Pendulous  Test  Mass _ 

Hammer  Sensitivity 


Model  No _ l0?fl07  _  Serial  No _ 7607  Sens.  9.91  mv/g 


^  •  05 _ lb  (__lZ^_gram)  including  accelerometer 


CONFIGURATION 

Tip 

PLRSTIC/VINYL 

PLRSTIC/VINYL 

Extender 

NONE 

STEEL 

SCALING  FACTOR 

Ib/g 

1.05 

0.99 

(SENSITIVITY  RATIO) 

(SEE  NOTE) 

Accel/Force 

IN/ms7) 

0.17 

0.45 

mV/lb 

5.48 

9.99 

HAMMER 

SENSITIVITY 

(mV/NI 

2.  13 

2.25 

NOTES: 

1.  The  sensitivity  ratio  <Sa/$f)  is  the  scaling  factor  for  convening  structural  transfer  measurements  into  engineeru 

units.  Divide  results  by  this  ratio. 

2.  Each  specific  hammer  configuration  has  a  different  sensitivity.  The  difference  is  a  constant  percentage,  which 
depends  on  the  mass  of  the  cap  and  tip  assembly  relative  to  the  total  mass  of  the  head.  Calibrating  the  specific  ham 
structure  being  used  automatically  compensates  for  mass  effects. 


Effective  mass 


0.265 


with  302A07  attached  and  vinyl  capped  plastic  tip 


PCB  PIEZOTRONICS.  INC  34?5  WALDEN  AVENUE  DEPEW.  NEW  YOR*  14043  ?495  'EitPHONt  716  684  0001  TWX  7’0  2 


Figure  62.  Force  Hammer  with  Calibration  Certificate 
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C&llnifeiraftndlli  (C®r(tfilfn<cai(l© 


Model  No.  336CB4 _ 

Serial  No.  ie;z6 

PO  No. _  Customer _ 

Calibration  traceable  lo  NIST  thru  Project  No. _ 922/25 11B1-93 


Per  ISA  RP37  2 


ICP1  ACCELEROMETER 

with  built -in  electronic  b 

Calibration  procedure  ts  In  compliance  with 
MiL-StU  45662A  and  traceable  to  NISI 


CALIBRATION  DATA  KEY  SPECIFICATIONS 


Voltage  Sensitivity 

99.2 

mV/g 

Range 

50 

Transverse  Sensitivity 

<,5.0 

% 

Resolution 

0.00! 

Resonant  Frequency 

17 

KHz 

Temp.  Range 

0/+150 

Time  Constant  (NOB. 

)  0.5 

s 

Output  Bias  Level 

12.8 

V 

±9 

g 

°F 


METRIC  CONVERSIONS 

ms  9  »  0  102  g 

“C-5/9X  (*F  32) 


Frequency  Hz 

10 

15 

30 

50 

100 

300 

500 

1000 

2000 

Amplitude  Deviation  % 

1.4 

1.3 

l.G 

.7 

0.0 

-.7 

-.7 

0.0 

4.7 

V 

FREQUENCY  RESPONSE 


Figure  63.  Accelerometer  SN#  1 0226  with  Calibration  Sheet 


Model  No  336C84 
Serial  No  10847 


CffiflnlbraiftMjl?  (Cteirftnffkaift® 

Per  ISARP37.2 

ICP‘  ACCELEROMETER 

with  built -in  elect/ omci 


* 


PO  No. _  Customer 

Calibration  traceable  to  NIST  thru  Proiect  No.  822/25 1 181-93 
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Figure  64.  Accelerometer  SN#  10847  with  Calibration  Sheet 
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Figure  65.  Accelerometer  SN#  1 0866  with  Calibration  Sheet 
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Figure  66.  Accelerometer  SN#  10877  with  Calibration  Sheet 
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Figure  67.  Accelerometer  SN#  10860  with  Calibration  Sheet 


Model  No.  336CB4 _ 

Serial  No.  _  10060 _ 

PO  No. _  Customer _ 

Calibration  traceable  to  NIST  thru  Project  No. 


if'CsrtnlfkM® 


Per  ISA  RP37  2 


022/251 101-93 


ICP*  ACCELEROMETER 

with  buill-in  electfomcs 

Calibralion  procedure  is  m  compliance  with 
MU  STD  45662A  and  traceable  lo  Nis  I 


CALIBRATION  DATA 

KEY  SPECIFICATIONS 

Voltage  Sensitivity 

96.5 

mV/g 

Range 

50 

±9 

Transverse  Sensitivity 

i5.0 

% 

Resolution 

0.001 

9 

Resonant  Frequency 

17 

kHz 

Temp.  Range 

0/+150 

®F 

Time  Constant  (NOH. ) 

0.5 

s 

Output  Bias  Level 

9.4 

V 

ME  I RIC  CONVERSIONS 

ms  *  -  0  102 g 

•C  -  S/9  x  (“F  32) 


m 


Frequency  Hz 

10 

15 

30 

50 

100 

300 

500 

1000 

2000 

Amplitude  Deviation  % 

2.1 

1.8 

2.0 

.7 

0.0 

-.7 

-.8 

-.6 

2.4 

FREQUENCY  RESPONSE 


+  3dB 


Amplitude  g 

Deviation 

*  3  dB 

1 

Id 

Frequency  in  Hertz 

00 

Piezolronics,  Inc.  3425  Walden  Avenue  Depew.  NY  14043-2495  USA 

/I6484-OOOI 


Figure  68.  Accelerometer  SN#  1 0868  with  Calibration  Sheet 
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Figure  69.  Accelerometer  SN#  10851  with  Calibration  Sheet 
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Figure  70.  Accelerometer  SN#  1 1 798  with  Calibration  Sheet 
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Figure  71 .  Accelerometer  SN#  1 0874  with  Calibration  Sheet 
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Figure  72.  Accelerometer  SN#  10867  with  Calibration  Sheet 
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Figure  73.  Accelerometer  SN#  1 0858  with  Calibration  Sheet 
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Figure  74.  Accelerometer  SN#  1 0856  with  Calibration  Sheet 
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Figure  75.  Accelerometer  SN#  10857  with  Calibration  Sheet 
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Figure  76.  Accelerometer  SN#  1 0854  with  Calibration  Sheet 
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Figure  77.  Accelerometer  SN#  1 0859  with  Calibration  Sheet 
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Figure  78.  PCB  Shaker  Specifications  Sheet.  Source:  [16] 
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APPENDIX  C.  SUPPLEMENTARY  TABLES 


Table  8.  FEM-Generated  Natural  Frequency  Comparison 


Mode 

Theoretical 

Natural 

Freq  (Hz) 

HIFI  FEM 
Natural  Freq 
(Hz) 

Theory-HIFI 

FEM 

Difference 

(%) 

Simple  FEM 
Natural  Freq 
(Hz) 

HIFI-Simple 

FEM 

Difference 

(%) 

Transition  FEM 
Natural  Freq 
(Hz) 

HIFI- 

Transition 

FEM 

Difference  (%) 

RIGID 

0.00 

0.00 

- 

0.00 

- 

0.00 

- 

RIGID 

0.00 

0.00 

- 

0.00 

- 

0.00 

- 

1 

14.51 

14.51 

0.0001 

14.51 

0.0008 

14.51 

0.0006 

2 

40.00 

40.00 

0.0000 

40.00 

0.0054 

40.00 

0.0044 

3 

78.41 

78.41 

0.0000 

78.42 

0.0201 

78.42 

0.0173 

4 

129.61 

129.61 

0.0000 

129.69 

0.0597 

129.67 

0.0489 

5 

193.62 

193.62 

0.0000 

193.86 

0.1283 

193.81 

0.1005 

6 

270.42 

270.42 

0.0000 

271.01 

0.2178 

270.98 

0.2076 

7 

360.03 

360.03 

0.0000 

361.41 

0.3821 

361 .29 

0.3504 

8 

462.44 

462.44 

0.0000 

465.59 

0.6817 

464.78 

0.5064 

9 

577.65 

577.65 

0.0000 

583.42 

0.9993 

582.25 

0.7962 

10 

705.66 

705.66 

0.0000 

714.57 

1.2631 

712.91 

1 .0278 

11 

846.47 

846.47 

0.0000 

861 .33 

1 .7552 

858.11 

1 .3744 

12 

1000.08 

1000.08 

0.0000 

1024.34 

2.4256 

1022.66 

2.2570 

13 

1166.50 

1166.50 

0.0000 

1191.95 

2.1817 

1188.16 

1 .8568 

14 

1345.71 

1345.71 

0.0000 

1474.10 

9.5406 

1445.81 

7.4384 

95 


Table  9.  Natural  Frequency  Shift  for  Various  Damage  Scenarios 
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Table  10.  Orthogonal  Projection  Sensitivity  Library  Selections 


ABCs  Available:  No  Pins,  1  Pin 

Modes  Available:  1-14 

Composite  [S] 

Row 

ABC  Configuration 

ABC  Location 
(Node) 

Mode 

1 

1  Pin 

12 

8 

2 

1  Pin 

12 

6 

3 

1  Pin 

12 

2 

4 

1  Pin 

9 

1 

5 

1  Pin 

4 

2 

6 

1  Pin 

4 

6 

7 

1  Pin 

13 

9 

8 

1  Pin 

9 

4 

9 

1  Pin 

13 

3 

10 

1  Pin 

6 

1 

11 

1  Pin 

9 

5 

12 

1  Pin 

4 

11 

13 

1  Pin 

13 

5 

14 

1  Pin 

8 

3 

ABCs  Available:  No  Pins,  1  Pin,  2  Pins 

Modes  Available:  1-14 

Composite  [S] 

Row 

ABC  Configuration 

ABC  Location 
(Node) 

Mode 

1 

2  Pins 

6,  12 

3 

2 

2  Pins 

3,4 

4 

3 

2  Pins 

7,8 

2 

4 

2  Pins 

12,  13 

9 

5 

2  Pins 

9,  10 

2 

6 

2  Pins 

5,6 

2 

7 

2  Pins 

4,5 

11 

8 

2  Pins 

10,  11 

1 

9 

2  Pins 

12,  13 

4 

10 

2  Pins 

4,5 

3 

11 

2  Pins 

3,4 

7 

12 

2  Pins 

4,  10 

5 

13 

2  Pins 

6,7 

1 

14 

2  Pins 

1,5 

6 

Clean,  1  Pin,  2  Pins,  3 

ABCs  Available: 

Pins 

Modes  Available: 

1-14 

Composite  [S] 

Row 

ABC  Configuration 

ABC  Location 
(Node) 

Mode 

3  Pins 
3  Pins 
3  Pins 
3  Pins 
3  Pins 
3  Pins 
3  Pins 
3  Pins 
3  Pins 
3  Pins 
3  Pins 
3  Pins 
3  Pins 
3  Pins 


10,  11,  12 

6,  7,8 

2,  3,4 
1,  7,8 

4,  5,  11 
2,  10,  11 
8,  9,  10 
10,  11,  12 
10,  13,  15 
4,  12,  15 

3,  6,  8 
8,  11,  13 


3,  4,5 
3,  6,  12 


ABCs  Available: 
Modes  Available: 


Composite  [S] 
Row 


Clean,  1  Pin,  1  Point  Mass,  1  Spring 
1-14 


ABC  Configuration  (mA4*Pa) 


_ 1  Pin 

_ 1  Spring  (0.014349) 

_ 1  Spring  (0.0028698) 

_ 1  Spring  (0.0028698) 

_ 1  Spring  (0.0028698) 

_ 1  Spring  (0.011479) 

_ 1  Spring  (0.0086094) 

_ 1  Spring  (0.0028698) 

_ 1  Pin 

_ 1  Spring  (0.0028698) 

_ 1  Spring  (0.0057396) 

_ 1  Spring  (0.0028698) 

_ 1  Spring  (0.011479) 

1  Spring  (0.0028698) 


ABC  Location 
(Node) 


Mode 


9 
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ABCs  Available: 

Modes  Available: 

No  Pins,  1  Pin 

4-14 

Composite  [S] 

Row 

ABC  Configuration 

ABC  Location 
(Node) 

Mode 

1 

1  Pin 

14 

6 

2 

1  Pin 

6 

9 

3 

1  Pin 

7 

4 

4 

1  Pin 

4 

6 

5 

1  Pin 

9 

5 

6 

1  Pin 

10 

5 

7 

1  Pin 

12 

4 

8 

1  Pin 

5 

6 

9 

1  Pin 

9 

10 

10 

1  Pin 

11 

10 

11 

1  Pin 

14 

7 

12 

1  Pin 

3 

9 

13 

1  Pin 

11 

6 

14 

1  Pin 

11 

8 

ABCs  Available: 

No  Pins,  1  Pin 

Modes  Available: 

4-8 

Composite  [S] 

Row 

ABC  Configuration 

ABC  Location 
(Node) 

Mode 

1 

1  Pin 

2 

5 

2 

1  Pin 

7 

4 

3 

1  Pin 

9 

4 

4 

1  Pin 

12 

6 

5 

1  Pin 

4 

8 

6 

1  Pin 

13 

6 

7 

1  Pin 

2 

6 

8 

1  Pin 

2 

7 

9 

1  Pin 

13 

5 

10 

1  Pin 

5 

8 

11 

1  Pin 

1 

4 

12 

1  Pin 

11 

8 

13 

1  Pin 

9 

5 

14 

1  Pin 

14 

6 

ABCs  Available:  No  Pins,  1  Pin,  2  Pins 

Modes  Available:  4-8 

Composite  [S] 

Row 

ABC  Configuration 

ABC  Location 
(Node) 

Mode 

1 

2  Pins 

8,9 

5 

2 

2  Pins 

12,  13 

4 

3 

2  Pins 

8,  11 

6 

4 

2  Pins 

2,3 

6 

5 

2  Pins 

5,8 

6 

6 

2  Pins 

13,  15 

8 

7 

2  Pins 

4,5 

4 

8 

2  Pins 

7,8 

4 

9 

2  Pins 

9,  12 

5 

10 

2  Pins 

9,  10 

7 

11 

2  Pins 

1,7 

5 

12 

2  Pins 

1,  10 

8 

13 

2  Pins 

4,  12 

4 

14 

2  Pins 

1,3 

4 

102 


Table  1 1 .  Initial  F EM -Experimental  Eigenvalue  Percent  Differences 


Resolution  (Hz) 

0.3125 

0.3125 

0.3125 

0.08 

0.02 

0.3125 

0.08 

ABC 

1  Pin 

1  Pin 

1  Pin 

1  Pin 

1  Pin 

2  Pins 

2  Pins 

Modes 

1-14 

4-14 

4-8 

4-8 

4-8 

4-8 

4-8 

Eigenvector  Row  1 

3.4637 

2.9473 

2.6898 

2.8634 

2.9555 

2.3973 

2.0436 

Eigenvector  Row  2 

2.6867 

4.8052 

2.0781 

1.9283 

2.2504 

8.5608 

7.3151 

Eigenvector  Row  3 

-0.7587 

2.0781 

2.0669 

2.0799 

2.4567 

2.4638 

1.9281 

Eigenvector  Row  4 

3.1215 

2.6975 

2.6097 

2.5103 

3.0618 

3.3268 

Eigenvector  Row  5 

1.7494 

2.4109 

3.2100 

3.2178 

3.3033 

1.4703 

1.2834 

Eigenvector  Row  6 

3.1224 

2.2428 

2.8386 

2.8172 

2.9656 

2.7247 

4.5891 

Eigenvector  Row  7 

5.7367 

2.4394 

3.0714 

2.9541 

3.0365 

0.2300 

1.0581 

Eigenvector  Row  8 

2.0632 

2.3877 

2.8435 

2.8428 

3.1403 

0.4634 

1.5534 

Eigenvector  Row  9 

4.3355 

4.8020 

2.4715 

2.5591 

2.5908 

2.2149 

2.2450 

Eigenvector  Row  10 

-1.5346 

5.2326 

3.9911 

3.9735 

3.9072 

2.4571 

2.4083 

Eigenvector  Row  11 

2.4109 

2.8098 

2.0361 

2.3619 

2.2059 

1.9038 

1.8402 

Eigenvector  Row  12 

6.0172 

4.6801 

3.4372 

3.3107 

3.4674 

2.9499 

3.3735 

Eigenvector  Row  13 

2.4691 

2.5715 

2.4109 

2.2858 

2.4475 

1.2182 

1.4020 

Eigenvector  Row  14 

4.5403 

3.4372 

2.9473 

2.9634 

2.6871 

2.1961 

2.4763 
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Table  12.  Updated  FEM  Comparison 


12 

974.70 

1027. 19| 

5.39| 

1011. 66| 

3.79| 

|96  '066 

1.67| 

993.19| 

1.90| 

1004. 48| 

3.06| 

993. 62| 

1.94| 

00 

CD 

o 

o 

o 

T — 1 

2.70| 

992. 6l| 

1.84 

IT 

826.40 

866.95 

4.91 

LD 

°o 

ro 

lo 

00 

3.32 

837.49 

1.34 

■3- 

o 

CO 

rvj 

00 

ro 

00 

1.44 

848.59 

2.69 

ro 

LO 

00 

ro 

00 

1.47 

844.78 

\llT 

00 

ID 

l< 

ro 

00 

1.36 

o 

T— 1 

690.26 

719.89 

4.29 

709.01 

\lLT 

691.29 

0.15 

694.80 

0.66 

704.65 

00 

O 

rsi 

695.52 

0.76 

703.12 

1.86 

694.22 

0.57 

CD 

564.78 

587.17 

3.96 

578.29 

2.39 

567.69 

0.52 

569.03 

LO 

o 

573.84 

1.60 

569.32 

o 

00 

o 

572.14 

1.30 

LO 

LO 

00 

CD 

LO 

0.67 

00 

452.94 

468.68 

3.48 

461.60 

1.91 

453.77 

0.18 

454.18 

0.27 

456.92 

00 

00 

o 

454.02 

0.24 

455.96 

0.67 

453.64 

0.15 

352.46 

364.09 

3.30 

358.59 

1.74 

353.54 

0.31 

353.15 

0.19 

355.04 

ro 

o 

353.20 

0.21 

355.25 

0.79 

353.14 

0.19 

in 

264.96 

273.06 

ID 

q 

ro 

268.93 

1.50 

265.21 

60'0 

265.10 

LO 

o 

o 

266.02 

o 

o 

265.20 

CT) 

O 

o 

266.76 

00 

ID 

o 

265.23 

01 '0 

LD 

189.84 

195.31 

00 

00 

rsi 

192.36 

1.33 

190.03 

0.10 

189.85 

o 

o 

o 

189.63 

0.11 

189.66 

0.09 

189.76 

o 

o 

189.602 

0.13 

127.04 

130.68 

2.86 

128.70 

1.31 

127.70 

0.52 

127.58 

0.43 

127.00 

O 

o 

127.52 

0.37 

127.32 

|eto 

127.44 

0.32 

m 

76.78 

79.04 

2.94 

77.84 

1.38 

75.06 

2.24 

76.25 

o 

o 

77.49 

[E6T) 

76.42 

0.47 

78.19 

1.84 

76.11 

00 

o 

CN 

39.22 

40.32 

2.79 

39.71 

1.24 

39.01 

0.52 

38.92 

0.76 

40.21 

2.52 

38.92 

0.76 

40.16 

2.41 

38.90 

18  0 

14.26 

14.63 

2.57 

14.41 

1.02 

14.221 

0.28 

14.24 

0.14 

14.23 

0.20 

14.23 

0.21 

14.24 

0.12 

14.20j 

0.42 

Mode  No. 

Test  Freqs  (Hz) 

Initial  FEM  (Hz) 

|Test-l niti al  FEM  Delta  (%) 

Adjusted  FEM  (Hz) 

|Test-Adjusted  FEM  Delta  (%)  ] 

Upated  FEM  A  (Hz)  | 

|Test-Updated  FEM  A  Delta  (%)  | 

Updated  FEM  B  (Hz) 

|Test-Updated  FEM  B  Delta  (%)  | 

Updated  FEM  C  (Hz)  | 

|Test-Updated  FEM  C  Delta  (%)  | 

Updated  FEM  D  (Hz) 

|Test-Updated  FEM  D  Delta  (%)  | 

Updated  FEM  E  (Hz)  | 

|Test-Updated  FEM  E  Delta  (%)  | 

Updated  FEM  F  (Hz)  | 

Test-Updated  FEM  F  Delta  (%) 
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Table  13.  Averages  of  Eigenvalue  Absolute  Percent  Differences 


Modes 

1-3 

4-8 

9-12 

1-12 

Initial  FEM 

2.769 

3.116 

4.637 

3.536 

Adjusted  FEM 

1.216 

1.558 

3.055 

1.971 

Updated  FEM  A 

1.014 

0.241 

0.919 

0.660 

Updated  FEM  B 

0.529 

0.190 

1.186 

0.607 

Updated  FEM  C 

1.214 

0.432 

2.357 

1.269 

Updated  FEM  D 

0.482 

0.201 

1.243 

0.619 

Updated  FEM  E 

1.455 

0.481 

2.022 

1.238 

Updated  FEM  F 

0.698 

0.179 

1.111 

0.619 
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